The Interaction of Radiation and Matter: 

Quantum Theory 

I. Quantum Mechanics of harmonic Oscillators:i 

To prepare the ground for the quantization of the electromagnetic field, let us revisit the 
classical treatment of a simple harmonic oscillator with one degree of freedom. We 
start, with the assertion that the Hamiltonian of such an oscillator must have the form 

2m'- ■' 

where p and q are the cannonically conjugate variables. Consequently, the accordant 
Hamilton equations yield a pair of coupled, inhomogeneous, first-order equations ~ viz. 

dq d!H 

-=^=/>/m [I-2a] 

dp d9{ 2 

-^=- — =-m(0'q. [I-2b] 

Combining these two equations, we see that formal classical analysis of the harmonic 
oscillator leads to a second-order homogeneous equation of motion 



dt 

with the evident complete solution 



d q 2 2 

— =-« 1 [I-2c] 



(0) 

q{t)= qiO) cosco t+ — — sin CO t 

mto [1.3] 

-m CO ^(O) sinco?-i- p (O) cosco t 



William homselY s Radiation and Noise in Quantum Electronics, McGraw-Hill (1964) has an excellent set of 
pre-quantum optics reviews. 
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Equations [ I-2a ] and [ I-2b ] are decoupled by transforming to a new set of dynamic 
variables ~ viz.'^ 

V2a) ^ ' ^ 

[1-4] 

y 2a) 

q= , \a + a^l 

so that , [ 1-5 ] 



• r t 1 



In terms of these new variables the Hamiltonian becomes 

^}^ =03 a [1-6] 

Thus, Equations [ I-2a ] and [ I- 2b ] are transformed to the decoupled normal mode form 
~ i.e. 

da . >v . da d9-[ 

dt dt da ^ ^ 

da^ . + . da^ d!H + 

= J CO a I = — - — = -co a 

dt dt oa 

With this background, let us now turn to the quantum treatment of the harmonic oscillator. 
From the general theory discussed earlier, we make use of the classical analogy and 
postulate that the dynamic variables q and p become the Hermitian operators q and p 
which satisfy the commutation relations 



^ In classical theory, the a and variables are frequently referred to as normal mode amphtudes. 

R. Victor Jones, April 11, 2000 
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] = 
[q,p ]=ih 



[1-8] 



The Schrodinger representation of the energy eigenvectors — i.e. u^[q) = (^l^^) — satisfies 
the following representations of the energy eigenvalue equation 



[q\f^\E„)=E/,q\E„) 

Zm 



1 

2m 



^.2 " , 2 ,^2 2 

-n — 2+"^ ^ 1 
dq 



[I-9a] 
[I-9b] 

[I-9c] 



The solution to this equation is familiar as 



/mO) 



m(0 7I 



[I-lOa] 



where {^^m (o/h qj is one of the Hermite polynomials discussed in the set of lecture 

notes entitled On Classical Electromagnetic Fields and the corresponding eigenvalues are 
given by 

E^=h(i)^ = h(i){n + l/2) . [I-lOb] 

An alternative treatment of the harmonic oscillator provides a powerful set of new tools. 
Recall the transformation law between state vectors in the Schrodinger and Heisenberg 
pictures — viz. 



R. Victor Jones, April 11, 2000 
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= T {t,0) |v|/^(0)) = exp j I V|/^(0)) 

which means that the operators transform according to the rules 

p„{t) = ^ t(^,0)p,T(r,0) 

The Heisenberg equations of motion for and p h 

dqjj__d9{jj__ 
dt dpu 

dp f, 2 
dt dq„ ^" 

By analogy with Equation [ 1-3 ], we see that the time evolution 

qfi{t) = ^ \t,0) qs 'T (?,0) = cosa)?+-^ sinto? 

Pu{i) = ^ \t,Q>) ps CT {t,Q>)=-m(A qs sm(At+ p ^ cosco? 
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[I-U] 



[1-12] 



[I-13a] 
[I-13b] 



[1-14] 



is consistent with Equations [ I-13a ] and [ I-13b ]. By analogy with Equation [ 1-4 ] we 
introduce two critically important (non-Hermitian) operators ~ i.e. 



a = \ ^ \(£> q + i p /m] 
a^= j ^ \(>} q - i p / m] 



[1-15] 



or 
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2m to 



a +a 



4 



with the commutation relation 



Thus 



a, a 



2 L' 

t 



a -a 



= 1 



a'^a + 1/2 



[fAf+V2] 



[1-16] 



[1-17] 



[1-18] 



t 

where = a is an important (as yet uninterpreted) Hermitian operator (observable). 
The Heisenberg equations of motion become 



da^ _ 1 

dt i h 

dal _ 1 

dt i h 



[1-19] 



with the solutions 



a^(r) = T \t,Q) T (?,0) = exp(-/a) r) 
fl J(?) = T \t,Q) fl/ T (r,0) = fl/ exp(+j to 

Consider the following equivalent eigenvalues problems ~ viz.^ 

Oi \E')=E\E') and 9^ \n')=n\n) 



[1-20] 



[1-21] 



First note that 



3 It follows from Equation [ 1-18 ] that E' = ^tO \n + 1/2] 
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[fA^ ,a]=!J\C a— a 9^ = a^a a —a a^a = 



t 

a , a 



a =—a 



r 



a , a 



t 



or 



Ot a^=a\9i +1) 



[ I-22a ] 



[ I-22b ] 



Thus, we see that 

91 {a = a {91 -\)\n') = {n-\)a\n') = {n'-\)\n-\) 
9t ^'^\n'^=a^ {9i, +l)\n') = {n'+l) a^\n')^{n'+l)\n'+l) 



[1-23] 



That is, if \n') is an eigenket of the operator 9^ with eigenvalue n' , then a \n'^f is also an 
eigenket, but with the eigenvalue n' - \ . Similarly, if \n') is an eigenket of the operator 

fA^ with eigenvalue n' , then a is an eigenvector with eigenvalue {n +l) . Hence, 
the operators a and a ^ "generate" new eigenkets which have associated eigenvalues 
differing by one from eigenvalues associated with the original eigenkets. It follows that 
each eigenvalue of 9(^ differs from all other eigenvalues by an integer! 



Normalization of these eigenkets proceeds as follows: 
Suppose a \n)= c^\n-\) 



[1-24] 



Then 



n{n\n) = [n\hC \n) = ^ fl |n) = c„(n| a ' \n-l) 
= c„(n-l|a \n)* =\c„f{n-l\n-l) 



ti 



[1-25] 
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but we require {n\n) = {n-l\n-l)=l so that c^ = 'Jn. Therefore 



9i \n)= n \n) 

a |0)= 

a I n) = \ n — X) 



[1-26] 



which leads to the powerful state vector generating rule ~ viz. 



|0> 



[1-27] 



There is another interest connect between the two alternative treatments of a harmonic 
oscillator: we may use Equations [ 1-26 ] and [ 1-27 ] to obtain the coordinate 
representatives of the harmonic oscillator states. Using Equation [ 1-15 ] we see that 



h d 

())q+ 

m dq 



<?|0>=0 



[1-28] 



Therefore 



(?|0)=Co exp 



? 

2h 



[1-29] 



and with the equivalence 



t 



m 



2/zcO 



fi d 

(Hq 

m dq 



[1-30] 



R. Victor Jones, April 11, 2000 
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Equation [ 1-27 } yields 




i_r 



(0 q- 



nn 

m aq \ 



[1-31] 



which is the generator of the eigenf unction displayed in Equation [ I-lOa ]. 



The connection between bosons and harmonic oscillators; 

There are two distinct interpretations of the energy formula = ^CO (n + 1) ; 



a. The formula describes the energy eigenvalue which corresponds to a given 
eigenstate of a particular harmonic oscillator. 

b. The formula describes the total energy of a system of n identical noninteracting 
quanta (particles), each of which is in the same dynamic state. The energy of the 
particle state is ^co. 

Again we turn to Dirac"^ who rigorously articulates the equivalence of these two view 

points. To quote his conclusion " This means that the dynamic system 

consisting of an assembly of similar bosons is equivalent to the dynamic 
system consisting of a set of oscillators - the two systems are just the 
same system looked at from two different points of view. There is one 
oscillator associated with each independent boson state. We have here 
one of the most fundamental results of quantum mechanics, which 
enables a unification of the wave and corpuscular theories of light to be 

effected Thus a set of harmonic oscillators is equivalent to an 

assembly of bosons in stationary states with no interactions between 



In Sections 60 and 61 of P. A. M. Dirac's The Principles of Quantum Mechanics (Revised fourth edition), Oxford 
University Press (1967). 
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them. If an oscillator of the set is in its n'th quantum state, there are n' 
bosons in the associated boson state. " 



R. Victor Jones, April 11, 2000 
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II. CANONICAL Quantization of Electrodynamics: 

With the foregoing preparation, we are now in a position to apply the classical analogy or 
canonical quantization program to achieve the second quantization of the 
electromagnetic field.^ As our starting point and for reference, we, once again, set forth the 
vacuum or microscopic Maxwell's equations in the time domain: 

VxE(f,f) = -|^B(r,?) [n-la] 

V xB {r,t) = |io J (r,r) +e^^\L^ ^ E (r,r) [ ii-ib ] 

at 

V E(r,?)=p(r,?)/8o [H-lc] 

V-B(r,t)=0 [IMd] 

The canonical formulation of classical electrodynamics ( Jeans' Theorem) is most 
conveniently achieved in terms of the (magnetic) vector potential in the time domain ~ viz. 

B{i,t) = 'VxA{i,t) [n-2a] 

E(f,?)=-|^A(f,?)- V(p(?,?) [II-2b] 

so that 

V[ V • A(f, t)] - V' A(i, ^) + 7 ^ r) + 1^ (p(f, t)=[ij(i,t) [ ll-3a ] 



In common usage, the process of treating the cordinates q. and as quantized variables is called first 

quantization. Second quantization is the process of quantizing fields — say, A{r,t) — which have an 
infinite number of dequees of freedom. 
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-80 V-|^A(?,?)-EoV'(p(f,?)=p(?,r)w [II-3b] 

In QED (Quantum Electrodynamics) it is convenient and traditional to make use of the 
Coulomb gauge -- i.e. V- A(f,?)= — so that 

V' A (f, - — A (?,?) =-[i, J^(i, t) [ n-4a ] 

c ot 

V' (p (r, t) = -p (r, t)/e^ [ 11.4b ] 

where 3Jr,t)=J{r,t)-JJr,t) = J{r,t)-e^^^(p{r,t) is the so called transverse 

at 

current density. Since A(f , t) is completely determined by the transverse current density 
in the Coulomb gauge, electromagnetic problems become in a sense separable ~ i.e. 

The transverse field problem: 

V-ET(r,r)= 

VxEfi,t) = -[i,^il{i,t) [n-5a] 
of 

Vx'H^[r,t)= J^[r,t) + ^j-E^[r,t) 

The longitudinal field problem: 

V-Ej^r,t)=p{r,t)/e, 
jJ^r,t)=-e,-E^{r,t) 
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We turn now explicitly to a treatment of the free electromagnetic field -- formally 
the case of ij^,t)= wherein 

13^- 

V'a(?,?)--2^A(?,?) = [II-6] 
We look for solutions in the form 

^ V s 

Substituting into Equation [ 11-6 ] we obtain 



2^ a, ? U3 r ^ /A +C.C. =0 



[11-8] 



Thus we may apply separation of variables techniques and the original problem is 
divided in to two new and distinct problems ~ viz. solutions of the following set of 
equations: 

2 

V'ui^) + '^ni^) = Q [II-9] 



ait)+(ii\QLlt)=Q> [11-10] 

where the (0 ^'s are separation constants. The spatial equations allow us to treat the 
boundary value problem of the cavity or defining field space in whatever detail that might 
seem appropriate in a particular case.^ But the bottom line is that the boundary condition 



The U,|rj 's do not form a trae complete set of solutions since no longitudinal vector field can be expanded in 
terms of such divergentless functions. 
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taken together with the characteristics of the normal modes or eigenfunctions determine the 
to s's which may, in turn, be identified as the eigenfrequencies of the normal modes. 



Thus, we may now write 



E 



and 



^ V^O s 

HtI? , r) = — V X A (f,?) = 2^{a V xii + c.c.}. 

M-o ^ vM-o s 



[II-U] 



[11-12] 



The crucial step required in establishing Jeans' Theorem is the expansion of the 
instantaneous value of the stored electromagnetic energy in terms of the cavity modes ~ 



VIZ. 



8o|E(f,?)f+|lo|H(f,?)f 



dV 



[II- 13a] 



which, in light of Equations [ II-l 1 ] and [ 11-12 ], becomes 
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To proceed we need the value of the integral I I I v x u • V x u * j • Using a well 



cavity 

known vector identity J we obtain 




[Vxu3- VxuV] jy= J J U3x(vxuV)-c/5+ J J J u 3 • Vx(vxuV) c/V [II-14a] 



By boundary value arguments we may easily show that the first term on the RHS of this 
equation vanishes and by using an even more famiUar (famous) vector identity^ for a 
divergenceles field, we obtain 




[Vxii3-VxiiV]^V=-JJJ u^-V^u:-dV='^ JJJ u^-K-dV [«-14b] 
Therefore, Equation [ Il-lSb ] becomes 

\^(0)= ^ s« CLXt) a At) III u 3 • uV dV [ 11-15 ] 



From Equation [ 11- 14b ] we may also write 



7 Namely.that V-(XxVxY)= (v xX) • (v X y) - X -(v X V X y) 
^ Namely, that 

Xx(VxY)] =(VxX)-(VxY)-X-Vx(VxY) 

= (vxx)-(vxy)-x-v(v- y)+x-v^ Y 
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and, again, by boundary value arguments we may easily show that the resultant surface imtegral on the left ^ 
Therefore, we may take 

JJJ [a..fi:.]jy=5., [11-17] 

so that Equation [ 11-15 ] becomes 

{W{t))=-J^0,\a,a: [11-18] 

s 

which when compared to Equation [ 1-6 ]^ is effectively the content of Jean's 
Theorem with 



J — a^^a= \(Oq+ip/m\ 

V 2 ^ \2(£> L ^' J 

^ a; ^ a' = 1^ [(£) q- i p/m] 



[11-19] 



9 That is !}{ =(0 a 
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To accomplish the canonical quantization program, field variables are expressed as field 
operators by making the identification 



V2 "■^i^'"*'^'''^"'"'^'' 



which leads to the following set of field operators: 



[ 11-20 ] 



Hxl? , t)= cTI—^ \a Vxu 4r) + a Vxu' 



The Plane Wave Expansion of the Electromagnetic Field Hamiltonian: 

To be definite, we may write an explicit plane wave representation for the field as 



[11-21] 



A(?,?)=LL j \ e,, L ,exp[/(k^ ?-(0 ^ f)] + a^t e^p [-/(k , ■ ?-co , t) 

s a=l y ^ s ^0 



[ II-22a ] 
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orio 



A {i,t)=ljlj^J [^soit] ^s.(r) + cil{t) nl(i)^ [ II-22b ] 



where 



Usa(f) = ;^e3<, exp[jk3-?] and a;„(f)= e^^ exp[-j ■ f] [II-23a] 

^»sa(?) = ^Jsa(0) exp[-/a),r] and al{t) = al{0) exp[+i(ii , t] [II-23b] 

Of course, ^^5= |kj=(0 ^/c in all of these expansions. Further the electric field 
expansion is given by 



= • j".. exp[i k, f- i (0. (] - fl t(() exp[i k. f- i 00 . t]| 



[n-24a] 



where £ ^ = I and the magnetic field expansion by 



H(f . <)= it [t.x S„(-r)] - «.t(<) [k. X S;,(f)]| [ II-24b ] 



10 



This expansion is often written as ^(r, — E (r,?) + E (Tj^) where 



E^H?'^)= 'XLcsa ^s^^sa exp[/ k, • f- /co , t] and E*-^(f,?)=-iXj;e3„ Es^^sIW exp[/k. 
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In Ught of Equations [ 11-18 ] and [ 11-21 ], the Hamiltonian of the radiation field is 

2 r . . ^ 2 

I I £ 

{'} ;^=-^ /^=-oo 

where k^^j = — [ x + /^y + /^z] [ 11-26 ] 

^{/} = c I ^{/} I 

The electromagnetic momentum (Poynting vector divided by c^ ) is given classical by 

JJJ [ll-27a] 

cavity 

and in terms of the second quantization operators it becomes 

^ = ^LL^ "^{Z} ^^J}c + «{/}a ^^J}a| [ ll-27b ] 

{/} <j=i 
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Thus, the Fock or number states are eigenstates of both the energy and the momentum 
of the field. 
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III. REPRESENTATIONS OF PHOTON STATES 



1. Fock or "Number" States: .11 



As we have seen, the Fock or number states 



[III-l] 



are complete set eigenstates of an important group of commuting observables — viz- 

Reprise of Characteristics and Properties of Fock States: 

a. The expectation value of the number operator and the fractional 
uncertainty associated with a single Fock state: 

{n\9i \n)= n [ in-2a ] 

An = [ "uncertainty" ] = ^\(^n\9l ' \n)-{n\9l = [ ni-2b ] 

b. Expectation value of the fields associated with a single mode: 

For one mode Equations [ n-24a ] and [ n-24b ] reduce to 

E(f , t) = ie 2;|a exp[/ k f- /co - a ^ exp k f - /o) ? 



[ in-3a ] 



H(f , t) = i £ [k X e |a (?) exp [/ k ■ r- / (0 ? - a \t) exp [-/ k ■ r- / (O ? | [ III-3b ] 



1^ In what follows, for simplicity we drop the k subscripts on the operators and state vectors with the obvious 
meaning that ||n - }^ ^ | ' ^ k ^ ^ ' 

R. Victor Jones, May 2, 20W 
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where £ = 




(n|E|n)= 
(n|H|n)= 



[ III-4a ] 



AE= /{(n|E.E|n)-(n|E|n)^}= ^(n + i^^ 





c. Phase of field associated with single mode: 

To obtain something analogous to the classical theory we would like to separate the 
creation and destruction operators (and, thus, the electric and magnetic field 
operators) into a product of amplitude and phase operators. Following Susskind 
and Glogower,i2 we define a phase operator, 4> such that 



Defined in this way, the basic properties of the phase operator may be evaluated 
from known properties of the creation, destruction and number operators. 
Inverting, we obtain 




a'^ = exp(-j<D)(fAf + l) 



[111-5] 



Susskind, L. and Glogower, J., Physics, 1, 49 (1964) 



R. Victor Jones, May 2, 2(2© 
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/ iinr A-y^ [ni-6] 

exp(-i4))=a' (fAL + lJ 

and since a =0\C + l,it follows that 

exp(/ <!)) exp(-j ^) = l [ ni-7 ] 

but only in this order! Operating on number states with the phase operators, 
we obtain from Equation [ 1-26 ] 

exp(j<D)|n)= {fJ\C +1)'^' a\n)={fJ\C +l)"^^ (n)^^ |n-l)= 

i- . i . [ni-8] 

exF(-j4))|n)= a^{9t +l)"^Mn)= a ' (n+l)"^^ |n> = |n + l) 

Consequently, the only nonvanishing matrix elements of the phase operator 
are 

(n -l|exp(/ \n) = 1 

[ ni-9 ] 

{n + l|exp(-j 4>)|n) = 1 

The phase operators defined by Equation [ III-36 ] do have the felicitous or 
classically analogous property of revealing magnitude independent 
information, but unfortunately they are nonHermitian operators — i.e. 

{n - 1| exp(/ <I>) \n) (n |exp(/ <I>) \n - 1)* 

~ and, hence, cannot represent observables. However, they may be paired 
into operators that are observables ~ viz. 



R. Victor Jones, May 2, 2®$ 
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cos<I> = i {exp(/ <I>) + exp(-/ <!>)} 
sin4> = ^ |exp(i 4>) - exp(-j 4>)} 



[ III-IO ] 



which have the following nonvanishing matrix elements: 

(n -l|cos4>|n) = (n|cos4>|n -1) = ^ 

(n-l|sin4>|n) = -(n|sin45|n-l) = — 

2 / 



[III-U] 



These nearly commuting operators may be adopted as the quantum mechanical 
operators which represent (as we will demonstrate anon) the observable phase 
properties of the electromagnetic field. 

For the Fock state: 

(n|cos4>|n) = (n| sin<I> |n>=0 [in-12a] 
Acos4> = Asin4>= J{<n|cos' 4>|n) cos4)|n)'} = JJ^ [ III- 12b ] 

Acos4) Asin4) = [III-12c] 

c. The coordinate or Schrodinger representation of state: 

Recall from Equations [ I-lOa ] and [1-31] that 



1^ Also, it may be easily established that the matrix elements of their commutator are given by 



(n|[cos4),sin4)]|n') = ^5„„, 5„o 



R. Victor Jones, May 2, 2(2® 
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Therefore, the probability P[q) of eigenvalues q for a given Fock state \n) is give 
by 

Df \ 1/ I M 1 Tj2 I /"JfO 1 r OT(0 2I rTTT 1/11 

d. Approximate "localization" of a photon: " 

Of course a plane wave is distributed or "de-localized" in both time and space. 
Defining the "wave function for a photon" is a task fraught with danger, 15 but the 
simpler task of defining a wave function approximately localized at a given instant 
is relatively straight forward ~ viz. 

' k- k ' 

2. Photon States of Well-defined Phase : 

Consider the state defined by 

|cp) = lim (5+1)"^ y exp[m (p]|n> [ III-15 ] 



See Section 10.4.2 in Leonard Mandel and Emil Wolf, Optical Coherence and Quantum Optics, Cambridge 
Press (1995), ISBN 0-521-4171 12. 

1^ See Section 1.5.4 in Marian O. Scully and M. Suhail Zubairy, Quantum Optics, Cambridge Press (1997), 
ISBN 0-521-43458. 

R. Victor Jones, May 2, 2(2® 
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Clearly, (cp |(p) = 1 given the orthonormal properties of the number states. Essential 
question: Is this state an eigenstate of the phase operators? To answer the question we 
need to consider the following potential eigenvalue equation: 

cos4>|(p) = - lim(5+l) ^exp[j n (p] exp[j4>] |n)+^exp[m (p] exp[-j 4>] |n) [ [III-16a] 
Using Equations [ III- 10 ] and [ III- 10 ], we obtain 



^ |(p) = - lim (s +1) yexpfm (pi \n - 1) + ^ exp[jn(p]|n +!)[ 

= - lim (5 +1)"^' j exp(/(p)^exp[/v (p]|v) + exp[-/(p] ^exp[iv (p] |v/ ' 



^ ..... ^j,., [III-16b] 

= cos(p |(p) 

+ - lim (s + 1)"^' |exp[/ 5 (p] |5 + 1) - exp[/ {s + 1) (pi |5> - exp[-/ (p] |0)| 

so that the state |(p) fails to be a strict eigenket of cos^> by terms that diminish faster 
than (5 + 1) as 5 . Similarly, we can see that diagonal matrix elements of cos^> 
and sin4> are given by 

((p|cos<I>|(p) = cos(p |l- lim (5 + 1)~'| => cos(p [ III-lTa 1 

((p |sin 4> |(p) = sin(p ] 1 - lim (5 + 1)~' [ sin(p [ ni-17b 1 



R. Victor Jones, May 2, 2(24 
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Reprise of Characteristics and Properties of Phase States: 

a. The expectation value of the number operator and the fractional 
uncertainty associated with a state of well-defined phase: 

sis +1) 



= lim — 



fractional 
uncertainty 



{<cp|^A^^|cp)-<cp|^A^ |cp)'} 



i 



iMmis+lY'Tn^- 



11=0 



lim (5+1) 



«=o 



lim (5+1) '52" 



«=o 



I lim !"-f2^^ + i')- — 



lim - 



73 



[III- 18a] 



[III- 18b] 



b. Expectation value of the fields associated with a single mode: 

From Equation [ III-3a ] 



((p|E|(p) =-2 I esin(kr-(0 cp) lim (^-l- 1) '^(n-l-l) 

diverges as Js for large s ! 



[III-19] 



R. Victor Jones, May 2, 2(2)9 
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c. Phase of field associated with single mode: 



((p| cos4>|(p) = coscp 
((p| sin4> |(p) =sin(p 



[ III-20a ] 



Acos4> = Asin4>= 



<(p|cos'4>|(p)-((p|cos4>|(p)'} =0 [III-20b] 



d. Probability of photon number: 

Finally, we may easily deduce the probability of finding n photons (i.e. the photon 
statistics) in a particular state of well defined phase — viz. 



We see that there is a equal, but small probability of any number: this agrees with the 
intuition that the magnitude of the field is completely undetermined if the phase is 
precisely known! 

3. Coherent Photon States: i6 

It would, indeed, be useful to have eigenstates of the destruction operator (electric 
or magnetic field) - viz. 



^„=Nrf-;mi(^+i) 



-1 



[ III-50 ] 



k/ 



[III-51 ] 



Reprise of Characteristics and Properties of Coherent States: 



a. The Fock state representation of the coherent state: 



The coherent state is a Harvard invention! See R. J. Glauber, Phys. Rev. 131, 2766 (1963). 
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Since, a ^\n)=yfn+l \n+l) and a = 9{^ +1 , then {n\ a = Jn+l {n+l\ and we 
are able to write a representative of the sought state in the number state basis ~ viz. 

{n\a \a)= •Jn+l{n+l\a)= a{n\a) [III-52a] 
a , , a" 



or (n|a)=-^<n-l|a)=-p<0|a) 
Using the expansion of the identity operator, the eigenket becomes 

|a> = I|n)Ha) = <0|a)i:^|n). 



To normalize the eigenket write 

(a|a) = <a|0)<0|a)i;^^ 



(a|0)(0|a>exp[|a|']: 



r 1 



Finally, we see that 



so that ( a| 0) = (0| a) = exp ^- - |a| 

|a)=exi^--|a| j2^-7^ K 



[ ni-52b ] 



[ III-53 ] 



[ III-54 ] 



[ III-55 ] 



is a normalize representation of the eigenkets of the destruction operator. 
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b. The expectation value of the number operator and the fractional 
uncertainty associated with a coherent state: 

(a^|a) = bc|' [ni-56a] 



[uncertaintyj (a^ |a) laNl "^^^l^nl 



||2n 
a 



[ ni-56b ] 



=iar 



Thus, we see that the fractional uncertainty diminishes with mean photon number! 

c. Expectation value of the electric field associated with a single mode: 

From Equation [ ni-3a ] 



(a|E|a) = -2 /-^ e|a| sin(k-r-a) [in-57a] 
y 28q V 

where a= |a| exp(/'d) . 

AE=J{(a|E.E|a)-(a|E|a)^}= [ni-57b] 



l'^ Similarly AH= j ^ ^ for the coherent state, so that AE AH =ch(i)/2V . 

cii,\2e,V 
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d. Probability of photon number: 

From the representation of the coherent state given in Equation [ III-55 ] we may 
easily deduce the probability of finding n photons (the photon statistics) in a 
particular coherent state is given by a Poisson distribution characterized by the 
mean value n = |a| . ~ viz. 



2n 



[ ni-58 ] 



Sample Poisson Distributions - Coherent State Photon Statistics 



0-4t r\. n = kr = i 




0.2t 



J 


i = |(x.p = i 






n = |o:.p = 10 






\ -„. 

1 n = |(x.| =^0 














n = |K.p = l 






-1 



40. 



SO. 



120. 
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e. Phase of field associated with single mode: 



(a|cos4>|a) = ^ exp 



1 

2 exp 



1 



- lafJE 



a a" + a a 



= lalcosiJ exp 



^1 



2« 



[ III-59a] 



Unfortunately, it is not possible to evaluate this summation analytically. However, 
Carruthersi^ has given an asymptotic expansion which is valid for a large mean 
number of photons ~ viz. 



(a|cos<I>|a) =cosi5 



1 



1 



7+. 



|ar»l 



[ III-59b] 



f. Coherent states as a basis: 

As we will see presently, the coherent states are very useful in describing the 
quantized electromagnetic field, but, alas, there is a complication — the coherent 
states are not truly orthogonal! From Equation [ III-6 ] we see that 



(P|a) = expr-i|af-i|pr1 l^t^ 
= exp[-l|af-i|P|Var^ 



[ ni-60 ] 



so that 



18 Carruthers, P. and Nieto, M. M., Phys. Rev. Lett. 14, 387 (1965) 



R. Victor Jones, May 2, 2(8)0 



The Interaction of Radiation and Matter: Quantum Theory 



Page 31 



(a|3)(P|a> = exp(-H'-|P|' + aP*+a*p) 

= exp(-(a-|3){a*-|3*))=exi(-|a-|3f) 



[ III-61 ] 



That is, the eigenkets are approximately orthogonal only when |a-P| is 
large! 

g. The "displacement operator:" 

There are a growing and significant set of applications where it is useful to express 
the coherent states directly in terms of the vacuum state |0). If we use the number 
state generating rule 



\n) = 



a 



t 



|0) 



~ i.e. Equation [ 1-27 ] ~ the coherent state may be written in the form 

( 1 2^ (oc^^^j r t 1 2^ 

|a)=exp--|cc|' £^ -^|0)=expaa'''--|cc|' |0) [III-62] 



If we make us of the Baker-Hausdorff theorem, we may easily show that 



The Baker-Hausdorff theorem or identity may be stated as 

exp{ + !B } = exp{ } exp{!B } exp{- ^ [ ,!S ]} 

when , ,!S ]] = [(S , , !B ]] = 0. For a proof, see, for example, Charles P. Slichter's Principles of 

Magnetic Resonance, Appendix A or William Louisell's Radiation and Noise in Quantum Electronics. 
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\a}= {a)\0)= exp(^a a - a* a j|0) 



[ ni-63 ] 



4. 

SO that ' (a) may be interpreted as a creation operator which generates a 

4. 

coherent State from the vacuum. (Its adjoint operator {d)= (-a) is a 

destruction operator which destroys a state). In some treatments ^ (cx) is 
described as the "displacement operator" (written D [a) )20 and the coherent states 
are called the "displaced states of the vacuum." 



To explore this point of view (and to give some meaning to the phase of the 
coherent state eigenvalue), we may express |a) in a two-dimensional, 

dimensionless "phase space" representation. To that end, following Equation 
[ 1-16 ], we write the dimensionless coordinate as 



e = 



^2ma)^^^ 



n 



q = exp [/ y] + a exp[-/ y] 



[ in-64a ] 



and the dimensionless momentum as 



% = 



p = exp[j(Y+7c/2)] + a exp[-i(Y+7c/2)] [ III-64b ] 



so that 



[e,7i]= 2i 



a , a 



= 2i 



[ in-64c ] 



We can (or rather you will) show that 2) ' (a) a 2? (a) = a + a and 25 ' (a) a ' ® (a) = a ' + a 

See Elements of Quantum Optics, Pierre Meystre and Murray Sargent III, Spinger-Verlag (1991), ISBN 0-387- 
54190-X. 
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and since these variables are canonical 22 

((AG)' )((A7t)' )>1 [ni-64d] 

Since 



t 1 

a ' = - (0- i Tl) exp[-/ y] 
a = - (0+/7r) exp[/ y] 



[ ni-65 ] 



the mode field (see Equation [II-24a]) b 

E[r, t)= ie'E^a exp[/k-r- /o)?]-fl ^ exp[-/ k-r+Zco ?]| [in-66a] 



becomes 



E(r,r)=-e £{71 cos(k r-co ?+Y) + e sin (k-r-co ?+y)} [III-66b] 

Since p has a coordinate space representation -ihd/d(J = -i[h(i>f2) ^d/dQ 
and ^ has a momentum representation ihd/d p = i [ft/ 2 0))^' d/dn, ^3 

aa^ - a* a =a^ -a^+ia^ + 

=-[a^ d/dQ+a. d/dn] [III-67a] 



22 



Of course, in general ((AA)' )((AB)' )>-^[ A,B] )|' where ( (A A)' )=( A' )-<A )' 



23 If this unfamiliar, see Equations [ 1-20 ] and [ 1-22 ] in the lecture notes entitled The Interaction of Radiation 
and Matter: Semiclassical Theory. 
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and 



A ^ (oc) = exp|^a a ^ - a* a j= exp[-(a^ d/Je+a, J/Jn)] ^ lll-67b ] 

Thus, fA^ (a) defines or generates a two-dimensional Taylor expansion when it 

acts on a function of 6 and 71. In particular, if we take the "phase space" 
representation of the ground or vacuum state |9jt) as the product of two Gaussians 

(see Equations [ I- 10a ] and [ 1-29 ]), then A ' (a)|0ji) represents a shift or 

displacement of this "phase space" representation ~ i.e. 

(eji|a) = (0jt|;l ^(«)|0)= «c(9-«J u^{j^-<x^ [ni-68] 
In light of Equation [ n-23b ], |a |a exp(- j co we can write 

(ejt|a(r))= Me(e-|a|cos((o r-t-c]))) Me(jt-|cc|cos((or-i-{|))) [III-69] 
where a= | a | exp(/ (|)) . 

h. The diagonal coherent-state representation of the density operator 
(Glauber-Sudarshan P-representation) : 



It may be easily established that 

71 



1 = - 11 IP) <PM' P= 11 IP) (PI d Re(P) d Im(P) [ ni-70 ] 



so that it seems quite reseasonable to look for a representation of the density matrix 
is the form 

P= JJ^(P)IP)<PM'P [ni-71] 

For a pure coherent state, P is clearly a two-dimensional delta function 
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Example 1 -- Coherent state 

^(13)= 5*'' ( p- a ) =5 ( Re (p) - Re (a) ) 5<" ( Im (p) - Im(a) ) [ III-72 ] 

In general, using Equation [ 111-60 ] - i.e. 

r 1 2 1 2 *1 

(P|a) = exp^- - |a| - - PI +ap J [ ni-60 ] 

we may find a simple procedure for finding the P-representation by writing 
(-a|p|a)=JJp(p)<-a|P)(P|a)J^P 

= exp(-|af) P) exp(-p|')jexp[ap*-p a*] / p 

Thus, (-a I p I a) exp |- |cx| j is the two-dimensional Fourier transform Of the 
function /'( P ) exp (-|Pf ) and we may write 

p( p ) = -L exp( Ipf) jj (-a| p |a) exp( |o(|') exp[-ap*+p a*] d'' a [ III-74 ] 

As a second example, consider a thermal radiation field described by a canonical 
ensemble 



exp(-i/f/^,r) 
Tr[exp(-^/fc,r)] 



[III-75 ] 



where ^ = h&\a'^ a + ■ Thus, 



P = I|1- 

n L 



exp 



exp 



n h(i) 



\n){n\ 



[III-76 ] 
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and 



{n)= Tr^p a a j= 2^|^exp 



\ 1 

1 



[III-77 ] 



so that 



„ (l + <n)) 



n) (n 



Thus, we can write (n | p | n) ^ 



and 



(! + (.))' 



in)" 



{-oc| P |«)=L , \' ^1 (-a|n)(n|a> 

exp(-^l') (-^f) r /„) T 



l + <n) 



exp 



(-Ml . J /r, , 1 11 



exp 



-^1 



V \"//J 



[III-78 ] 



[III-79 ] 



[III-80 ] 



Finally, we see that 

Example 2 -- Thermal radiation - a chaotic state 



( M1 



r 



Tl\\ + {n)) [ 
71 ^ ' 



1 



^1 



exp(-|3a* + a|3*) |3 



[III-81 ] 



As a third example, consider Fock or number state. From Equation [ III-55 ] we 
see that 
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exp(-|af) „ 
(-a|p |a) = <-a|/i)<«|a)= ^ ^ -|a| [III-82a] 



and 



P( P )= li, exp( IPI^) jj(-H')" exp[-ar+P a*] a 

/,,2\ [III-82b] 

so that 

Example 3 -- Pure Fock or number state 

exp( Bn ;a2« 

n\ 3(3 3|3 



i. The Glauber-Sudarshan-Klauder "optical equivalence" theorem: 

Suppose we have some "normally ordered" function 



f(n) 



a , a 



[III-83 ] 



The expectation value is given by 

/f<^^[a,a'^|=Tr 



(N) 



r 



a , a 



[III-84 ] 
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Using Equation [ III-7 1 ] we see that 



/f(N) 



a , a 



Tr 



« m 

= jj^(«)EE^«.(«|a'"a"»^'a 

n m 
n m 

or, finally, the "optical equivalence" theorem 

[M a, a^\=\\p [a] f^''\a, a) 



V )l 

j. The Uncertainty Relationship for {0,7t} 



Since 



a , a 



t' 



■ 1 we see from Equation [ III-64a ] that 



/Afl2\^/n2\_/n\2 



l^a^ a^^l exp[2/ i\+{a a) exp[-2 ij] + {a ^ a^+^a a 



[III-85a ] 



[III-85b ] 



-|a exp[2/ y]-(^J/^ exp[2/ y]- 2^a ^^(a) 



[ III-86 ] 



= (:Ae':)+l 



where (: A :) symboUizes the normally ordered expectation value of the operator 
A . From Equation [III-85b ] 

/:Ae' a,a'^^}^i= \\ P{a) ^Q\a,a)d' a [III-87] 



V 
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^: A ' j^a , a ''"'j :^ = jj P (a) [a a* exp {iy)+Aa exp(-/ y)] ' a [ III-88 ] 



7 



/ r 

If we choose y (and P( a) ) such that ^:A0^ 



,.tl:\<0. 



then (A0^)>1 and 



( A ^> l(squeezed states)! 
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IV. INTERACTION HAMILTONIAN -- COUPLING FIELDS AND CHARGES: 23 

To build a complete quantum picture of the interaction of matter and radiation our first and 
most critical task is to construct a reliable Lagrangian-Hamiltonian formulation of the 
problem. In this treatment, we will confine ourselves to a nonrelativistic view which, 
fortunately, is adequate for most circumstances. We start with a representation of the 
Lorentz force for a single charged particle ~ viz. 

f=^E + vxB} [IV-la] 
or in terms of the electromagnetic potentials 

f= ^- V(P-^ + vx(vxA) I [IV-lb] 

Let us now write the total time derivative of a component of the vector potential 



dt 

so that 



« v-V U+— ^ IV-2] 

dr^ dt dt ^ ' dt 



vx(VxA) = v(v- a)-(v-v)a = v(v-A)--^ + |^ [IV-3] 



Therefore, we may write the Lorentz force as 



Much of what follows draws heavily on material in Chapter 5 of Rodney Loudon's The Quantum Theory of 
Light (second edition), Oxford (1983). Marian O. Scully and M. Suhail Zubairy in Section 5.1 of Quantum 
Optics, Cambridge (1997) have a slightly different, but quite insightful treatment of the subject of the atom- 
field interaction Hamiltonian. 
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5 r 1 

We may also write ~ A = - — Up - Vg Ao 



so that Equation [ IV-4] becomes 



d d d [ I 1 J a r 1 2I 

= m — v„ = m 1 - Vr Vr 1= 1 - mv , 



[IV-5] 



which may, in turn, be written 



This last equation may now be compared to the Lagrangian equation of motion ~ i.e. 

^\^L(i,y)\-^L (i,y)= [IV-7a] 

where, in general, 

L (f , v) = [Kinetic Energy] - [Potential Energy] = {x:,y)-U (r, v) . [ lV-7b ] 
Therefore, we identify 

X(f,v)=^^mv'-^(p-v-A) j [IV-8] 

as the Lagrangian for a single charged particle. We may write 

(?, v) = -^X (?, v) v„+-^X {r,y) v„ [ IV-9a ] 

ar„ ' av„ 
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which, in light of Equation [ IV-8 ], becomes 



l-L (r,v) = (r,v) v„+-^i: (r,v) v„ 



[ IV-9b ] 



Therefore, 

V - — L (f , v)-X (f ,v) = A constant of motion [ IV-9c ] 

and, using the cannonically conjugate momenta associated Equation [ IV-8 ] ~ i.e. 

^a= (?>v)= mv„+^A„ -, [IV-10] 

ov 

we see that Equation [ rV-9c ] can be written 

- [^a - ? A J - ^ [ - ^ A J [ !P„ - ^ A J + qi?- 1 A„ - ^ A J ,,,,,,, 
m 2m m [IV-lla] 

= A constant of motion 

or finally 

1 r - - i2 

q(p+ — \'P-qA\ = A constant of motion . [IV-llb] 

We are now in a position to set forth the nonrelativistic Hamiltonian of a single 
charged particle ~ viz. 



y{ = q(p+ — \ ^ -qlY [IV-12] 
2m'- ■' 
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with the canonical conjugate variables given by 

^9-[ 1 

5(p 1 r 1 5 r i 

In principal, we are done, since we may now write the complete Hamiltonian for a many 
particle material system as 

^=L^, = LU,-9+^[2'>^a]' [IV-14] 



Unfortunately, this form of the Hamiltonian is not the most useful in optical physics since it 
is expressed in terms of the vector potential and needs to be evaluated at all charge 
positions. Most annoyingly, it does not yield an interaction term in the form used earlier - 
- i.e. Equation [ III-3a ] in the lecture set entitled The Interaction of Radiation and 
Matter: Semiclassical Theory. 

The Multipole Expansion of the Classical Hamiltonian: 

Before we manipulate the Hamiltonian further, we digress a bit to get a better fix on what 
we are really look for. From the most basic notions of electrostatics, the energy of 
interaction with an external transverse field should be expressible as the "energy of 
assembly" viz. 

\ ?=0 fj 

all C all C all C 

J E^(?) df=-1^^, J E^(f) df+1^^, J E^(f)-df [IV-15a] 



^■^ This is the work required to assemble the atom in the external field and of course ignore all inter-particle 
electrostatic interactions 
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which becomes for a neutral atom with charge centered at f = 

f/.t=£?, [IV-15b] 



Electric fields of interest vary only slightly over an atom so that we should be able to 
expand the external field in a rapidly converging series as follows 



E,(f) =E,(0) + (f- V) E,(f)^ + i(f. V)' E,(f)|_^^ + ... [IV-16] 



Substituting this expansion into Equation [ IV- 15b ] and integrating we obtain 



electrons 



u,u=-e L ir^yi+^,(rr^yj,(rrf + --]H^^^ [IV-17a] 

Formally, we may cast this expansion for the interaction energy in the form 

demons fexp(r. -Vj-l ^ , .. ^ 

In the continuum picture, the interaction energy should be expressible in the form 

fffEJf)-PmdV [IV-18] 



where P(r) is the polarization density. It may be shown quite easily that Equations 
[ IV-18 ] and [ IV- 17 ] are equivalent if the polarization density is expressed in the 
following 
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multipole expansion 

electrons 



p(f)=-e J2 {f-l/2!r(r- V)+l/3!r(r- vf +--j5(r-rj 



[IV- 19] 



With this background, we return to a consideration of the Hamiltonian in Equation 
[ IV- 14 ]. In what follows, we show that this Hamiltonian is transformed into a form 
which consistent with Equation [ IV- 17b ] if we make an appropriate gauge 
transformation of the fields. In general, the gauge transformation 



A.(r,t) =^ A(r,r)- V%(r,r) 
(p(r,f) ^ (p(r,?) + ^^ %(?,?) 



[ IV-20 ] 



where is an arbitrary scalar gauge function, leaves the electric and magnetic fields 

unchanged. Motivated by the discussion above, we choose the gauge function 

[ IV-21 ] 



ir,t)={l/e) A(F).p(F)dV' 



where P(r) is the polarization density in the form shown in Equation [ IV-19 ] and A(r) is 
the vector potential appearing in Equation [ IV- 14 ]. The impact of the transformation on 
the scalar potential term in Equation [ IV- 14 ] is easy to evaluate — viz. 



exp 



f V 



[ IV-22 ] 



Dealing with the transformation 
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A{r„t) ^ A(tO-(V^) jjjv,A(r).p(f)dV 



[ IV-23 ] 



is straightforward, but extremely tedious (and not very useful)! Substituting for P(f) 
from Equation [ IV- 19 ] and doing a lot of integrating by parts we could demonstrate that 



[ IV-24 ] 



Making use of these results, the complete Hamiltonian of a atom plus radiation field may be 
written 



2m 



L (w) J 



1 



[ IV-25 ] 



IL 



r.V 



J 



J 



In practical terms, the useful complete Hamiltonian is written 

= !}{ ^ + !H ^+[!H + ^ EQ-^ ^ MD-^ ^ Nl) 



[ IV-26a ] 



where 



= E "Pi "hrn- i e Y^a^v) +^-Ze (p(0) 



[ IV-26b ] 
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2JJJ KE/+|ioH']dy 

i 

^Ea=\eL v{ V] }=-V a- 



^Mo="-^L [I X 'P] .n{Q)='-^M . H(0) 



2m 



^ NL = 



8m 



L[^xH(0)]' 



[ IV-26C ] 
[ IV-26d ] 

[ IV-26e ] 

[IV-26f ] 

[IV-26g] 



If we take i" of the order of magnitude of the Bohr radius a-^, M of order h, and 
V ET(r)|__^ of order |k| |Et(0)| then we may estimate the three linear interaction terms^^ ~ 



VIZ. 



me 



^ EQ ~ Et(0) 



2) eh 
16 mc 



e h 



^Mz,-H(0)^ = E,(0) 

2m 2 mc 



Further we note that 



^ and ^ 



1 



ED 



ED 



4ne^hc 137 



25 We assume here that Hd} =h c k= e^jla^ = R^/2 
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V. PHOTON ABSORPTION AND EMISSION 

"POOR MAN'S" Second Quantization of Material System: 

In treating the complete quantum mechanical problem, it is useful to recast the material 
(atomic) Hamiltonian in terms of an appropriate set of creation and destruction 
operators. To that end we make the following definition 

^\x)= n(a^\x) [V-1] 
Using the ubiquitous identity operation ^|j^X-^l = 1 > we may write the material 

X 

Hamiltonian in second/ quantized form ~ viz. 

= i:\y){y\=Il.V) ^co. {x\y){y\ 

X 

= L^0). |x>(x| 



X y 



[V-2] 



In general, the operator |x) (y| = applied to any state \z) yields 

\^){yV)=KK\z)=\^)K 



[V-3] 



— i.e. the operator changes a state z to a state x if the state is y otherwise it produces zero. 
In other words, the operator destroys the state y and creates a state x. The second 
quantization viewpoint is particularly useful in treating the interaction of a two-level 
material system with the radiation field. This case, is most conveniently formulate in two- 
vector notation with the use of Pauli spin matrices — viz. 



and \by- 



{a\ = [l 0] and {b\ = [0 l] 



[ V-4a ] 



[ V-4b ] 
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\a){a\ 



"1" 


[0 


1]= 


"0 


r 











"o" 




0] = 


"0 


0' 


1 


[1 


1 






= o 







0] = 


"1 


0" 





[1 








"0" 


[0 


1]= 


"0 


0" 


1 





1 



Consequently, the atomic Hamiltonian may be written 



"1 0" 




"0 0" 













1 



1 
-1 



1 

2 V a 1 



and if we neglect the mean energy of the states 



2 



ab 



1 
-1 
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[ V-4c ] 



[ V-4d ] 



[V-5a] 



[V-5b] 
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and the electric dipole interaction Hamiltonian becomes 

= [I«)<«I+I^)H]M][I«)<«I+I^><^I]Et(o) 

= {|a) (a| [e + [e^D]|a)(a| } Et(0) 

= {e{a\'D \b) c^+e(b\'D \a)c~]-E^{0) 
= -{ po^ +p o~ }-Et(0) 



[V-6] 



From Equation [ II-24a ] in this lecture set we can write 



{;} <3=i ■> 



where •£{,} = 



'{/} . 



2e,V 



is the so called the electric field per photon and f, is the 



location of the center of the atom under consideration. Thus Equation [ V-6 ] may be 
written quite generally for a two level atom as 



^ED= {e{a\'D \b) (5^+e{b\'D |a>a"} 



X 



{01 *=i 



{ i ^{1} s{ t) exp[j k{,} rj^-i a^^ ^ t) exp[-/ k^,} • r, ] } [ V-8 ] 



where the coupling constant is given by 
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8u}s= ^ ■ e{.}, = e {a\ '3\b)- e^,j, [ V-9 ] 

Interaction of a Two-level Atom and a Single Mode Field -- Rabi Flopping; 

Let us first consider the interaction of a two-level system with a single photon state to 
make contact with the discussion in Section in of the lecture set entitled The 
Interaction of Radiation and Matter: Semiclassical Theory. Equation [ V-8 ] then 
reduces to 

^ED = i ^ {si + sl exi(/ k • r,] - al{t) exp[-j k • f^] } [ V-lOa ] 



where g^= h e2^p^-e-=e I (a| ^D | ^) • e- . If we neglect any inter- atomic 

W 2 ^ Cq 

interference effects by taking = 0, we may simplify Equation [ V-lOa ] to 

^ED = ^{ 8i +4 }{ iHi^)-^ } 

= n{g-^G^+gl<y-}{ia-^{t) + adj.} [V-lOb] 
Thus, we may then write the complete effective Hamiltonian of the composite system as 
IH: =!H, + !H,= Ui<Si^, o^ + ^a)-aV+^{ gj^o"+g^o-}{ / a-^) } [V-lla] 

In the rotating wave approximation this reduces to 

!H =9{, + !H,^Uiiii^ o' + h(ii-^a'a+h{ig-^ao^-i^^a' o-} [V-llb] 
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In the spirit of the discussion in the Section VII, Semiconductor Photonics in the 
lecture set entitled The Interaction of Radiation and Matter: Semiclassical Theory, 
this effective Hamiltonian may be adapted to provide a fully quantum mechanical 
treatment of optical interactions in semiconductors .^^ 

"Dressed" Atomic States: 

We know that the unperturbed Hamiltonian satisfies the following eigenvalue equations 

^K^\an)=ti^^ i(0,, + n(0-j|a n) 

r 1 1 ^^"''^ 

^,\bn)=r^-^(ii^, + n(i)-^'^\bn) 

- where 1'^ I'^/I^r) and |^")-|^)hk) -- and the electric dipole perturbation 
couples the states \a n) and |Z7(«+l)^. It is useful to resolve the complete 



2^ In Semiconductor Photonics we noted that Chow, Koch and Sargent in their Semiconductor-Laser Physics 
(Springer- Verlag - 1994) treat semiconductor problems in terms of the following semiclassical Hamiltonian for 
an inhomogeneous two-level system: 



where k' ^ k} and il' ^ k} are, respectively, electron and hole {creatiop destructio^ operators. M-k is 

the dipole matrix element between vertical states in the valence and conduction bands. In the fully quantal 
treatment the effective Hamiltonian becomes 

where ^- = >/2 |I- sin(fc,. z^)/fl 
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Hamiltonian into a sum of component Hamiltonians ^ =^^„ where the component 



9{ act only within the { |a n) , \b (« + !)) } coupled manifold of states and can be 
writterf'? 



'1 0" 








.0 1. 




"1 0" 


n\ 




+- 


.0 1. 


2 



n\ 5(0 2 / a 1 
2[-2 igla' -60) J 



5a) 



2 / g- Jn+1 1 



-2 ig*- -Jn+l 



-5a) 



[V-13] 



where 5a)= a)^ ~fi>i;. The second term in this equation is of the same form as the 
coupling matrix in Equation [ III-8c ] of the lecture set entitled The Interaction of 
Radiation and Matter: Semiclassical Theory where the semiclassical Rabi frequency 
Q.^ = p Eq/H is replaced by its quantum equivalent ~ viz. ^ i ^ Jn+l . DiagonaUzing 
this matruc we find the eigenvalues of the so call dressed atomic states -- viz.^^ 



2^ In particular, using the identity operator - i.e.,. / = T^, |t|/ \T|| 

^[\an){an\+\b{n + l)){b{n + l)\] {\an)(an\ + \b{n + l))(b{n + l)\] 
= \an)(an\ {an\9{ \an) + \anYb{n + \)\ {an\9< \b{n + \)) 

+ |&(n+l)Xan| [b{n+\)\j{ \an) + \b{n + \)Yb{n+\)\ {b{n+\)\^ \b{n+\)) 

2^ It should be noted that the rotation matrix 

^ rcos0„ -sin0„ 
sin9„ cos9„ 

diagonalizes the Hamiltonian in Equation [ V-13 ] through the transformation 9t i?/^ ^ Further, 9t 

relates \hedressed andbare probabihty amphtudes as 
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r v-141 

(see energy level diagram below) where ^„ - ^5(0 +4|^| {n+l) is the quantized 
field generalization of the Rabi flopping frequency. 



The dressed eigenstates are given by 

|ln)= sin6„ |a n)+cos0„ |& (n + l)) 
\2n)= cos9„ |a n)- sin6„ |& (n + l)) 



[ V-15a] 



where cosQ„= . ^ ^ [V-15b] 

JK-H +4|g,|(n+l) 

and sine„= ■ ^ ^ [V-15c] 



n 

= LK{t)\^n)+C,Xt)\2n)] 
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The time evolution of a state is directly represented in terms of these dressed states 
viz. 



|x|/(r))=exp(-/i?/' rM)|v(0)) 

cxj 2 

= L£exp(- iE^„ t/n]\s n^sn\^\f{0)} [ V-16a ] 

n=0 s=l 
oo 

= Lexp(-/[n + i]a)-,r){ exp(-i iQ^ t)\l n)Q„{0) + ex^i iQ^ t)\2 n)C,„{0)} 



n=0 



or more explicitly 



rc,„(r)1^rexp(i/n^) 1rQ„(0)1 

Cu{t)\ [ exp(-i/n^)J[Q„((|J 



[ V-16b ] 



Therefore, the flopping of the undressed states is given by 



ex^-iiQ^t)\ K^i)(0). 



9t 



rcos( i t) + i sin{ \ t) cos2e„ -/ sin{ I t) sin20„ 1 [ CjO) 1 

-/ sin( i t2« r) sin2e„ cos( i r) - / sin{ i t) cos2d„\[c,^^,){0) \ 



[V-17] 



_ "cos(ia«?)-/(5a)M«) sin{ia^) -i(2g^^l/a^) sin{ ia« If Q„(0) 1 
-/(2g^^/^/^2:) sin(in^) cos( + / (60)/^:) sin( ?)j LQ-^i)((|J 
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Perhaps the most reveahng application of the this result is for the case of a resonant 
coupled system ~ i.e. 6co= 0— which is prepared so that C'^„+i)(0) = 0. In this 

instance. Equation [ V-17 ] yields 



which clearly exhibits the simplest manifestation of spon^aneoMS emission — i.e. 
Rabi flopping in the absence of an applied field! ! ! 

Interaction of a Two-level Atom and a Multi Mode Field -- Spontaneous 
Emission; 

To broaden (make more realistic) our treatment of spontaneous emission we return to 
Equation [ V-8 ] to include the interaction with many modes with a two-level atom ~ 
viz. 



This equation may be easily generalized to encompass multi-level material systems.^^ 

If we are dealing with situation in which the locations of the atoms are uncorrelated 
we may, for simplicity, dispense with the exp[±j k^^j -r^ factors — i.e. we will 

neglect, for the present, any possible interference effects -- and write 
For a multilevel atomic system the effective interaction Hamiltonian in second quantized form becomes 



|C.„(Of=cos^[^-7^r]|C,„(0)f 



[V-18] 




{/} s=l 




CCP {/} 5=1 V 
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{1} s=\ 



i(^{t)s{t)-ic^l}kt)\ [V-19] 



Transforming to the Schrodinger picture and taking the unperturbed ground state 
|^) = |£»)|0)|0)---|0) as the zero energy reference point we may write the complete 

effective Hamiltonian as 



2 , ^ 



{f) s=l 



[ V-20a ] 



If we include only energy-conserving terms ~ i.e. in the rotating wave 
approximation ~ 



2 J. 



^" ]} ■ 



[ V-20b ] 



It is important to note that, in general, the interaction terms in this Hamiltonian include 
contributions from the coupling of the material system (atom) to any externally excited 
mode(s) (the incident electromagnetic field) and to all available electromagnetic cavity 
modes. For the present, we ignore the coupling to externally excited modes: we treat 
the external interaction later as a perturbation. Our goal at this point is to diagonalize 
the Hamiltonian of the complete unperturbed system which may be written 
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[to, 





. * 

* 

2 



[-igl 
= • (c5 + g) • a 





igl ■ 






• 


■ 


1 ' 







• 


\ \ 



CO, 



a 



[V-21] 



We are looking for the eigenstates and eigenvalues of this Hamiltonian which we write 
as 

!H\s)=hWs). [V-22] 

Since the square matrix is Hermitian, it is possible, in principle, to diagonalize it via 
a unitary transformation of the form 

U-(ro+g)-U"' = ^ [V-23] 

where A, is a diagonal matrix whose elements are the A'^+l eigenvalues of the 
dressed states of the coupled system. If we define the row vector = a^ • and 
the column vector jl = U • a then Equation [ V-22 ] becomes 

N+\ 

h~^!H = 2L^-{& + 'ij-2L = ^^ ■X-^ = Y,\\^\\^t [V-24] 

For consistency, (r||i.||i.j |5)=5^j 5^^ and, hence, and |I are, respectively, paired 

creation and destruction operators in the sense of the operators defined in Equation 
[ V-3 ] above ~ i.e. \s)=\i\ \g) and |^)=M., I^). From Equation [ V-23 ] we may 

write 
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which can be written explicitly as 





U,2 










U23 • 




U31 


U32 


U33 • 





U(Ar+i)i ^[N+l)2 ^{N+l)3 



1r (Oo ''^1 
1 1 ■ * 

2(^+1) 

-ig* 



U 



■igN 



I, 


Uu 




^1 Ui3 




U21 


K U22 


A,2 U23 


A,3 


U31 


A,3U32 


^3^33 



.^(iV+l)U(iv+l)l \N+il^{N+i)2 ^(Ar+iU(^+])3 



Expanding out the matrix product on the left we see that 



for elements in the first column and 



for elements not in the first column. Hence 
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[ V-25a ] 



[ V-25b ] 



[ V-26a ] 



[ V-26b ] 



[ V-26c ] 
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and all elements of the unitary matrix can be expressed in terms of the first-column 
elements as 



r 



u= 



U 

u 



21 



31 



L 



u 



A,2— cOj 



(w+^i 8i 



A,2— CO2 

i U31 g2 

^13-0)2 



{N+t)l Si 



^'Uu Sn 1 
^' U21 8n 

i U31 



(w+i)i Sv 



J 



[ V-27 ] 



Quite generally the columns of any unitary matrix satisfy an orthonormal condition 
viz. 



[ V-28a ] 



so the normalization of the first column gives 



[ V-28b ] 



and the orthogonality of the first column with any other column yields 



■= 



for each one of the cavity frequencies o) . 



Further we may multiply Equation [ V-28c ] by a product of factors 



[ V-29c ] 
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I P 

L;7^-n(co,-\)=0 [V-28d] 

It is obvious from this expression the co^.'s are roots of the left side of the equation if 
the to^.'s are replaced by CO . Thus 

n (w- K) = n (w-O) J [ V-28e ] 

, to— A, , 



Finally, we see that 

.2 n(co-coJ 
■-— [V-29] 



, 11(0)- 



We can make use of this expression to obtain the time varying polarization induced by 
an externally excited field. The Hamiltonian associated with this perturbation may be 
written 

^ED = fi{8e., +8L <^~}{iae.,- ^ ^l} [ V-30 ] 

Since a^= jl^ • U and a= U~' • jl we see that 

(J" = &!u,i and o" = £u;V,=Lu>, [V-31] 

t t t 

and from Equation [ V-6 ] we may write 

'D = (a\'D\b) +{b\'D\a)G- 

= !{<«! 'D \b) U, iil + {b\ S I a) u: n J [ ^-32 ] 
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In Ught of Equation [ V-29 ], the standardized form for the frequency dependent 
susceptibility (see Equation [ VA-12 ] becomes 



X{oy)=^{a\^)\b)\ \4— - + 



[ V-33 ] 



By eliminating the U's from Equations [ V-24a ] and [ V-24b ] we see that 



[V-34] 



Again 



^s-wo-Ll 8t r(^,-«r)"'jn(?i,-coj= 



so that we can write 



«-«o-Lkr(^.-«r)' 



[V-35] 



n(to-toJ=n(a)-Xj [V-36] 



Therefore 



e 



he 



«-«o-LUt f K-«)"' CO+(Oo-£|gj'(cO,+Co) I 

t t J 

where the sum ^| I (^r gives an explicit, non-phenomenological accounting 

t 

of interactions with the cavity modes and hence of spontaneous emission! 



i \ [ V-37 ] 
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Evaluation of Spontaneous Emission Rate: 

Recall the discussion of phenomenologically defined damping in semiclassical models 
of the dielectric response function in the lecture set entitled The Interaction of 
Radiation and Matter: Semiclassical Theory. Recollect, in particular, Equation 
[ III-19c ] in those notes. In reconciling that discussion with the content of Equation 
[ V-37 ] above, we see that the summation ^| g, \ (co, -(o) replaces the simple 



damping parameter y . Our task here is evaluate this integral which we write as 

I |2 f I P 

/ Y (co) = limV, — = lim — p,, ((O.) Jo), . 

^ ' , (0,-(0- /£ *' (0,-(0- Z8 ' 



[ V-38 ] 



where p„,(co,) Jtt), is number of cavity modes with frequency between (o^ and 
CO, + Jol),-- Following arguments best explicated long ago by Heitler,^'' it may be 
shown that 



W. Heitler, in Chapter II, Section 8 of The Quantum Theory of Radiation (3rd edition), Oxford Press (1954) 
uses contour integral arguments to shown that 

lim — ^- — =fPr - + i'Kh{x) 

e^O X- i£ X 

To quote W. H. Louisell's summary in Chapter 5 of Radiation and Noise in Quantum Electronics, 
". . .if /(x) is well behaved and has no poles at JC = a , then 



\l^dx = Tr \^dx + inf{x) 

•J X — a J X— n 



X- a -> X— a 

where 'Pr means the Cauchy principle part and is defined by 



J X- a Ji^M J X- a *' X- a 
provided the limit on the right side exists." 



J 
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r \ gf f 2 

to, —CO 





which we write as ij {(£>) = A{(£)) + i . This is an extremely important 

result! ! It shows that the interaction between the atom and the cavity modes leads to a 
frequency shift or correction in the atomic splitting to^^, 



A((o) = fPr -^^p,,((0,)rf(0, [V-40a] 
to, -co 

' 

and a spontaneous emission decay rate 

Y(co)=7iUj'p„,(co) . [V-40b] 

If we assume that the cavity modes defined for the blackbody calculation in Section IV 
of the lecture set entitled The Interaction of Radiation and Matter: Semiclassical 
Theory (see Equation [ IV-5 ] in those notes) are the appropriate modes, we know that 



pJco)=co7c' 
and from Equation [ V-9 ] we know that 

Treating the shift A(to) , the radiative correction to atomic energy level separation, is a 

very complex and much studied matter. The simple interpretation of Equation [ V-38a ] 
is problematic since the integrand is proportional to co^ at large co^ and, thus, the 

correction significantly diverges! ! 
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The divergence in ^to) was for many years an unresolved discrepancy between the 
quantum theory of radiation and observational spectroscopy. The difficulty was 
overcome by Bethe in 1947^1 using a technique known as mass renormaUzation. 
Bethe point out that the divergence can mainly be associated with the mass of the 
electron. It is found that the energy of a free electron has an infinite contribution arising 
from the interaction of the electron with the electromagnetic field. In other words, the 
apparent mass of the electron is shifted by an infinite amount from the mass of an 
electron which is not in interaction with the radiation field. However, the former mass 
is the one measured experimentally, since it is never possible to isolate an electron from 
the radiation field. Identification of the measured electron mass with the theoretical 
mass, after renormaUzation to take account of the energy of interaction with the 
radiation field, removes most of the divergence from ^co) 

Calculations for the hydrogen atom show that ^co) vanish unless one of the states in 
the transition is an S state. Even when it does not vanish, the renormalized A(co) is 
always very small compared with the excitation frequency (o^, and varies slowly with 
CO. For example, the magnitude of A(co) for the 2^5^ state of hydrogen is about 

10^ Hz, or roughly six orders of magnitude smaller than the n = 2 state excitation 

energy The existence of level shifts was first demonstrated by Lamb and 

Retherford in experiments on radiative transition between the 2 state of hydrogen 

and the unshifted 2 ^Py^ state. The splitting between these states is known as the Lamb 
shift. 32 



31 Bethe, H. A., Phys. Rev. 72, 339 (1947) 

32 From Chapter 8, Rodney Loudon, Quantum Theory of Light (1st edition), Oxford (1973) 
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However Equation [ V-38b ] is not complicated by divergences and, consequently, we 
easily obtain the famous Weisskopf-Wigner formuP^ for the spontaneous emission 
decay rate into the modes of a three-dimensional cavity 

^co) = n k\\H= ^^^^ \ia\ 'D \b). e,f 
^ 1 2(0^ 2 



33 V. Weisskopf and E. Wigner, Z. Phys., 63, 54 (1930). 
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APPENDIX: THE DIELECTRIC SUSCEPTIBILITY 
A General Dressed State Formulation 

Suppose that the complete Hamiltonian of a coupled system is parsed into two components 

i^ = i^„ + i^„. [VA-1] 

The component includes the Hamiltonians for the unperturbed material system, the 

free radiation field and interactions of the material system with available cavity modes. 
The component i?/"^ is the Hamiltonian for the interactions which couple the material 

system to externally excited modes. As the first step in finding a fully quantal 
expression for the dielectric susceptibility, let us expand the state vector in the Schrodinger 
picture in terms of, presumably, known eigenkets of ff{ ^ — viz. the dressed states of 

the unperturbed system ~ 

|¥(?)) = Lcir)exp(-/\f)|5> . [VA-2] 

s 

Following a now familiar track, we can use the Schrodinger equation of motion - i,e, 
to obtain 

" . [ VA-4 ] 

-/ hC:{t) = l^Clit) cx^+i (\- K) t] {q\lHJr) 



In turn, we obtain the following expansion for the time dependent expectation value of 
induced material system dipole moment: 
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(p(^)) = -(vWk^n¥(^)) 

= exp[ i K t] {r\ e 'D £c,(?) exp[-i t] \s) [ VA-5 ] 

r s 
r s 

Differentiating this expression with respect to time and using Equation [ VA-4 ] we obtain 
( P(^) )= -H)^ KlO Cit) exF[ / (\- g t\ <?|i^ Jr)| ^(r| S \s) 

+ Il[-Ic:(OC,(Oexp[-/(\-gr] 'A^ M)\l{r\'^ V) [VA-6a] 

r s 

Regrouping, we see that this expression can be written 

(^(^))= -Ilj^ Kl^) Clt) exp[ / [X^- X) t] [<?|i^ JrXrl 2) |.)- <r|^,J.X?l |r)]j 
- H.C:{t)Clt) X) exp[ X) t]e{r\S\s) 

r s 

Since 9-[ is proportional to !Z> , we see that, like magic, the first term vanishes! ! ! 
Hence, 

( P(^) ) = -H.<{t)Cit) i {X, - \) exp[/ [X^ - X) t] e{r\ ^D |.> . [ VA-6c ] 

r s 

Differentiating this expression with respect to time and, again, using Equation [ VA-4 ] we 
obtain 



[ VA-6b ] 
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( P(^) )= -IJllKitKit) - exp[ t] 

r s q ^ 

X [(X,-Xj(^|i^Jr)(r|^D|.) + {^,-Xj(^|2) IrXrKJ.)] . [VA-7] 
- IlClOCiO {K - KY exp[ / {X^ - y t] e{r\ 'D \s) 

r s 

Our task is to now to attempt an interpretation this very nasty expression. To that end, we 
make use of Equation [ V-33 ] to write Equation [ VA-6c ] as 

- - • [VA-8a] 

Xe{r\l^[{a\^\b}U,,^, + {b\^ |a) U>,} |.> 

t 

Using the properties of the and operators (v/z. |5')=|J,j|^) and |^)=|J,j5)), this 
expression reduces to 

{ P(0 ) = -ILc:{i)CXt) exp[ i{K - K) t] e {{a\ 6 \ b) U,, 5,^ +(&| 6 \a) 5, J [ VA-8b ] 

r s 

which may interpreted as a sum of a series of dipole moment components - viz. 

(p(^)) = Lp.(^)=-LK(^K,(?)exp[/\?] e{a\'D \b)\5,,+ c.c] [VA-8c] 

s s 

Given this interpretation, we return to Equation [ VA-7 ] and use Equation [ V-6 ] to obtain 

( m )= +ILLc:{tK{t) I [2K-K- K] X exp[ K) t] UO)-(q\ kXr| 'D \s) 

' ' " .[VA-9] 
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Again using Equation [ V-33 ] and the properties of the jx^ and jx^ operators, we see that 
{q\S\r){r\S \s)= (^^{(al U,,^I + <Z,| |a) Ut^^Jlr) 

X 

x<di:{<a| 2) 1^) U^, \i\ + {b\ S I a) U>4 |.) 

y 

= i:i:{<«l^""l^>'u.,U^,<^||lI|rXr||i;|.) 



+ 



la\S\b) 



U.,U^,<'?|^i:k)<r|^J.) 
^{b\'D\af\J:, U:,<^|nJrXr|nJ.) 
-f]<&|!P|a^'u:.U,,<^|nJrXr|n;|.) 



[VA-lOa] 



which reduces to 



(^|!Z)|rXr|2)|5>=(a|!Z)|&)(&|!Z)|a>U,,U:5, +(&|2)|aXa|2)|&)|U,J 5, 5, FVA-lOb] 



Substituting this expression and the expression in Equation [ V-33 ] into Equation [ VA-9 ] 
and, again, using the properties of the \\,] and operators it relatively straightforward to 

obtain 



( p(^) yi.^] vkt)= -'-iixk^k\ cimt) u^, u: 



q s 



xexp[i(>.^-\)?] ET(0)-<fl|!P IbXbl'D |fl).[VA-lla] 
+ ^ J^2^^ \c^{t)f I U,,f E,{0)-{b\ S \a){a\ S \b) 
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Thus, in a kind of rotating field approximation, we get a set of driven harmonic 
oscillator equations of the form 



Vit) + 'k]^lt) ==-E,(0)•|(a|^Z)|&) 



\^M-\^i'i Kf- [VA-llb] 



In the Weisskopf-Wigner approximation! ~ i.e. \C^{t)\ ~\ --we can easily solve these 

equations and sum their results to obtain a standardized form for the frequency 
dependent of the dressed dielectric susceptibility of the system 



«-4f N^^M' Et^I" 



N 



2 ,.^2| rl| 
1 



■ + ■ 



1 



A, -to A, +0) 



[ VA-12 ] 



1 V. Weisskopf and E. Wigner, Z. Phys., 63, 54 (1930). 
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VI. PHOTOELECTRIC DETECTION AND PHOTON COUNTING 

To introduce the subject let us quote Glauber in Photon Statistics^'^ 

"Photon counters function by absorbing photons from the field. What they count 
are, strictly speaking, not photons but atomic photoabsorption processes rendered 
individuall detectable by an amplification mechanism of some sort. In most of the 
devices used experimental the fundamental absorption process is the photoemission 
of an electron. The amplifying mechanism for the photoelectrons is usually a 
cascade multiplier. The precise way in which amplification is achieved need not 
concern us here however, since it plays no direct role in determining what the 
counter detects. 

"An ideal photon counter would be one which is quite small in size and has a 
sensitivity independent of photon frequencies (at least over the spectral of the 
incident field). Both requirements can be met rather well, in fact, by the simplest 
possible model of a counter, a single atom which we observe to see whether and 
when it emits a photoelectron . Since the atom is quite small compared with the 
wave length of visible light most of the transitions it can undergo may be treated by 
means of the electric dipole approximation, in this approximation the atom is 
coupled to the field through the interaction Hamiltonian 

^ =-eLqy-E(f,?) 

Y 

in which is the spatial coordinate of the yth electron of the atom relative to its 
nucleus which is located at r . 



From Fundamental Problems in Statistical Mechanics U (edited by E. G. D. Cohen), North-Holland 
Publishing (1968), pp. 140-187. 
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"The intensity of a field is proportional to the rate at which a counter records 
photons, or in the case of our single-atom counter to the probability per unit time of 
our observing a photoabsorption process. ..." 



Electrode 




With careful designed, high gain photomultipliers, it is possible to register the atomic 
ionization caused by individual, isolated photons. Following Glauber then we will take the 
photoelectric interaction Hamiltonian as 

^PE=-eY.\-'^[^^t) [VI-1] 

n 

where the r^'s are the relative spatial coordinates of the electrons bound to a nucleus 
located at R . First-order perturbation theory informs us that the transition amplitudes 
associated with the process of photoelectric absorption are proportial to the matrix elements 
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0=-^I</I?«-[e*''(R'^)+e*"'(r,?)]| [vi-2] 

n 

where / and i signify, respectively, the, final and initial states of the complete electron- 
photon system (See, for example. Equation [ 11-16 ] in the lecture set entitled The 
Interaction of Radiation and Matter: Semiclassical Theory. Recall from Equation 
[ II-24a ] that 

E^'(r,?)= J^E-q^q / a, exp[-/a)-?-h / q-R] 

' [ VI-3 ] 

e(-'(r,?)= L-E-qg-^l-ifl-;) exF[+/tO-?-/q-R] 
q 

If the system starts in the lowest electronic energy state and energy is conserved in the 
transistions, only destruction operators contribute in first-order and, hence, the rate of 
photoemissive transitions from one particular initial electronic state (i.e. the transistion 
probability of the detector atom absorbing a photon from the field at a position between 
times t and t+ At) is given by 

KE(^,t) = Y\ (/I f • E^'(r, t) I /) I' [ VI-4a ] 

Using Equation [ 1-6 ] ] in the lecture set entitled The Interaction of Radiation and 
Matter: Semiclassical Theory and the now familiar identity ^ li'X*! = 1, we see that 

= I( i I f • e"(r, ^) I / )( / I f ■ e<^'(r, ^) I / ) 

= {/|f.E<-'(R,r)E*^'(R,r).f|0 ^^'-^^^ 
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If we factorize the initial state as | / ) = | / = {/', = | A(i') ) | F[i") ) where | ^«') ) and 
I F[i") ) are, respectively, the initial electronic (atomic) and field (photon) states. 

w^R,t)= (A(0| (F(0| E*-'(R,r) E^'\R,t) |F(0) •f|A(0) [ VI-4c ] 

Summing over all initial states, we have the total probability of the detector atom absorbing 
a photon from the field at a position between times t and t+ dt. 

w..(R, t) = i;(A(r)| r |x;(F(r)| e<-'(r, t) E^^\R,t) \F{rf^ • r|A(r)) 

=i:(^(oif'-^"*'(^'^)-^^i^(o) [vi-5] 

where ^"®(r, r)= E<-'(r, r) E("'(R,r) 

i " 

Thus, the total rate of photoemission from a given atom depends on the expectation value 
of the operator^" ®(R,r) (or E^"'(r, E^^^'^R,?) ) and this rate, in turn, translates -- 

eventually ~ in the value of the observed photocurrent in the photomultiplier. 

(^"«(R,r)) = IPr(F(0)(F(0|E'-'(R,?)E*^'(R,^)|F(0) [VI-6a] 

However, p^^d^ ^^^^^(0) (^(^'01 write 

i " 

( ^""'(r,^) )= Tr[p,,,E<-'(R,?) e(^)(r,?)] [ VI-6b ] 
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We may then assert that this operator is (or is proportional to) the photon intensity 
operator which has an expectation value proportional to the observable light intensity. ^5 
and is a particular value of the more general function 



( G<"(r„ ? ■ R„ f,) ) = Tr[p,,, E<-'(r„ f) e'^'(r„?,) 



[VI-7] 



which is call the first-order correlation diadic (function). 



Hanbury Brown - Twist Effect — Correlation of Radiation from Incoherent Sources ^6 

If we extend the arguments associated with the introduction of Equation [ VI-4a ] to 
find the joint count probability of a double photexcitation 



,(r„ R„ = E| (/I ? ■ E<^'(r„ f,) E<^'(r„ ?J ■ ? I /) 



[ VI-8a ] 



Again using Equation [ T6 ] ] in the lecture set entitled The Interaction of Radiation and 

L|A/J = 1 
I / \ I ~ , we see that 



vv;,(R,,?-R,,f,)=(A(/)|r/(A,(/)|r;:(F(/)| 



E<-'(R„f,)E<-'(R„f,)E'^'(R„?,)E<^'(R„?,) :r, |A, (/)) rjA (/)) 



[ VI-8b ] 



This interpretation is consistent with Equation [ II-27b ] which gives the Poynting vector as 



c 



{'} o=l 



z 



{/} o=l 



Therefore, ( /(R,?) )= 2e„ C Tr[ p,.^, E<-'(r, ?) E*"'(r,?) 

36 R. Hanbury Brown and R. Q. Twiss, Proc. Roy. Soc, A 242, 300 (1957). 
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Summing over initial states 

w,,(R,j{,R,j,)=Y,{A,{i)\i;{A, (01 f;:/*''(Ri,?i; R,,?^) -r^ K (0) ?i K (0) [ ^I-Sc ] 

i 

where 

g^'\R,j,;R,j,)={F{i)\E^-\R,j)E^-\R^^^ [ VI-8c ] 

The Hanbury Brown - Twist Spectrometer 
I L 



Light 
Source 




Quite generally, we can follow Glauber an define the nth-order correlation function (tensor) 
as 
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Tr 



[VI-9] 



Quantum Interpretations of Fundamental Interference Experiments 3? 
Young Interference : 

To simplify things, suppose that we parse the radiation modes involved in the Young's 
experiment as follows: 

Mode 1: Light distribution corresponding to the emergence of photons from 

pinhole 1 . The field operators associated with this mode are i ' ^ i'} . 
Mode 2: Light distribution corresponding to the emergence of photons from 
pinhole 2. The field operators associated with this mode are 2 ' ^ 2 } . 







Mode 1 





Mode 2 



3' Much of this draws heavily on excellent discussions Chapter 6 of Rodney Loudon's The Quantum Theory of 
Light (2nd edition), Oxford University Press. 
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Suppose that the two pinholes provide the only means of escape for the photons. The 
number of photons n, which would be measured by a phototube directly behind ith 
pinhole is given by the expectation value of a/ a, .If the pinholes are of equal size, then 

a = ( fl. + a ,)/ Jl. 

/ ^ [ VI-10] 

(the [>/2 j factor insures that all operators satisfy the same commutation relationship). 
It is relative straight forward to show that 

^ [VI-11] 

[1 m i 7 r "I wi ifl r i m— 1 

a\ a. = a. \a\ +y^[«J 

where _/ = 1 or 2. Suppose that the incident state (with n photons) can be written 

l«> = [^/^^[«TlO> [VI-12] 

so that 

(n|fl;aj«) = l(0|[a]"a;ajaf|0) 

= ;^(o|«/[«]1«T«.|o)+^|<o|[a]""VT"|o) 

= -^<«-lk-l/=^ [VI-13] 
Thus for a n-photon, single mode incident beam, the intensity at Q is given by 

[VI-14a] 

+2 C* cos[ A: (5^ - Sj)] I «2 1 
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or 

{I (e . « - ^ { ^.r + iQf + 2 c: c, cos[k{s, - } [ ] 

Thus, the same intensity distribution as observed for a large number of photons n can be 
built up by doing n experiments with one photon! The interference is a one photon 
effect and does not depend in any way on the interaction of photons with 
each other. 

Hanbury Brown- Twiss Interference: 




■ 

u 

Experimental Configuration Transmission and Reflection Modes 



Following an argument identical to that used above, we see that the correlation between 
photomultiplier currents is proportional to 

/It \ / V / I 1 1 I V n{n-l) 

\«i «2 «2 «i/ = \'^i'^2/=\'^ki «2 «2«ih/= — — [VI-15] 
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Thus, the degree of second-order coherence is 









g%)={rHn,)l 


^"'^"^^-(./2)(n/2) 


(n-1) 
n 


[ VI-16 ] 
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Bunched light beam 



Random light beam 



Antibunched light beam 
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VIA. ON CLASSICAL OPTICAL INTERFERENCE, COHERENCE AND 
FLUCTUATIONS 

Light from a real physical source is never, strictly speaking, monochromatic nor does it 
emanate from a single point in space. Both the amplitude and the phase of the wave field 
generated by a real source undergoes irregular fluctuations. Within a time interval of the 
order of the putative coherence time the amplitude and phase of the radiation remain 
relatively stable and, thus, it may be treated as monochromatic. For time intervals long 
compared to the coherence time, the effects of fluctuations become manifest. In particular, 
interference phenomena, the hallmark of wave behavior, are observed only when the 
interfering beams traverse paths that differ in length by an amount small compared to the 
coherence length — i.e. the product of the coherence time and the velocity of light. In 
this first section we review the classical analysis and interpretation of interference effects in 
terms of the notions of partial coherence . 

COMPLEX Representation of Polychromatic Fields: 

Let us first carefully define the mathematical framework of the subject. We assume that the 
observable optical field may be represented by a Fourier integral transform i 




[VIA-1 ] 



Since of necessity E' [t) must be real 



For this assumption to be valid, it is sufficient that the function E{t) be square-integrable — i.e. 
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£(-0)) = £^(a)) [ VIA-2 ] 

and, consequently, Equation [ VIA-1 ] may be re-expressed in terms of an integral over 
positive frequencies ~ viz. 

E^'\t)= ^ {E{oi)cx^-im] + c.c]d(>i ■ [VIA-3] 



Further, if we write 

£(a)) = ^ a(a))exp[/(t)(a))] , [VIA-4] 

where a(a)) and M,^) are both real, then the observable optical field may be expressed in 
the form 

E^'-\t) = I a(o)) cos[/ [^((o) -« t]}d(xi [ VIA-5 ] 



which is a strict generalization of the complex or phasor representation of real 
monochromatic fields. To complete the generalization, we follow Dennis Gabor^ and 
introduce the so called complex analytic signal 

E{t)= J a((o) exp{/[(^((o)-(0?]}ja)= 2 j£'((o) exp[-/a)?]Ja) [VIA-6] 



which is to be associated with the real observable field ~ viz. V) ~ ^^[^(O] • As we 
shall see, this rather careful, or should we say pedantic, introduction of the analytic signal 



2 D. Gabor, J. Inst. Elec. Engrs. (London), 93, 2529 (1946). 
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greatly facilitates and validates the use of complex quantities in handling cycle-averaged 
properties ~ e.g. energy and intensity ~ of polychromatic optical fields. ^ In most 
applications of interest in optics, the spectral amplitudes will have appreciable values only 
within a spectral width Aco which is small compared to some mean frequency (to) . Thus, 
it often useful to express the analytic signal in terms of a temporal modulation of a mean 
monochromatic ~ viz. 

E{t) = a(Oexp [i [(^{t) - (co) t]} . [ viA-7 ] 

Using Equation [ VIA-6 ] we see that 

a ( Oexp [/([)(?)] = 2 J £ (co) exp [- / (co- (co)) f] Jco = J £ (co) exp [- / (at] d(xi[ VIA-8 ] 

^co) 

where £(c»>) = 2 £'(co+ (co)) ; by assumption, £(o)) will be appreciable only near co= (co) . 
The following list of transforms are included for reference: 

a {t) = ^[E''\t)]\[E%)]" = p{t)E*{t) = \E{t)\ [ VIA-9a ] 



^ The following identities involving improper functions are invaluable when treating analytic signals — see Beran 
and Parrent: 

^ I exi(/ x{y'- y))dx= d{y'- y) 
j- |exp[+/x(/-y)] dx= 6,(/-3;) = i[5(/-y)Tl-^j 
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(])(?)= (a))r+tan 




f E*{t)-E{t) l 
[' E\t)-E{t)\ 



[ VIA-9b ] 



E^''\t) = a(t) cos[/ [(^{t) - (to) t]} 
E^'\t) = a(t) sm{i [<^{t) - (to) t]} , 



[VIA-9d] 



[ VIA-9C ] 



Using Equations [ VIA-1 ] and [ VIA-6 ], we may easily establish a form of Parseval's 
formula 



We have thus far assumed that field is defined for all values of t. In practical terms, the 
field or, more likely, the observation of the field is defined only within some finite time 
interval -T <t <T . In all instances of any consequence, T is much larger than the 
physically significant times 2jt/Ato and 2rt/(to) so that we may feel some confident in 
taking the idealization T , However, stationary^ requires that the time averaged 
energy of the field, which is proportional to 





[VIA- 10] 








A random process is characterized as stationary when all of the moments of the random variables are 
independent of absolute time. 
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T 



Um — 



dt ' 



[ VIA-11 ] 



-T 



must be finite and, thus, a stationary field function is not square-integrablel The problem 
of analyzing such functions has been a much discussed issue in the mathematics literature.^ 
However, as applied physicists, we brush these mathematical niceties under the rug in what 
follows. We proceed by firmly asserting that optical fields of interest do, indeed, have 
Fourier transforms and chance the consequences of trucation errors. 

The cycled- averaged intensity or power spectral density is a crucial element in the analysis 
of coherence effects and, from Equation [ VIA-6 ], it is proportional 



Consideration of this expression leads us to a logical definition of the first-order 
temporal correlation function of the electric field as 




[VIA- 12] 




T 




[VIA- 13] 



-T 



So that Equation [ VIA- 12 ] becomes 



5 See, in particular, N. Wiener, Acta. Math.,55, 1 17 (1930). 
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|£((o)f=^-^j|(£(r+x)£*(r))exp[/a)x] Jx= j j r(x) exp[/ wx] Jx [VIA-14 



and the normalized power spectral density (distribution) is defined as 



r 



F{(o)= \E{(o)\ I |£(co)| Jco 

LI J 



r 



1r, 2\1"' 

, {E{t+x)E*{t)) exp[i (Ox] dx \ {\E{t)\ ) 
1% 1^ J J L\ /_ 



1 



— I y"'(x) exp[/o)x] dx 



[VIA- 15] 



where 



Y®(x). 



{E{t)Er{t)) 



[VIA- 16] 



This normalized correlation function is usually called the field's complex degree of 
first-order temporal coherence 

Young's Interference Experiment -- A Rudimentary Measurement of 
Temporal Coherence 

By considering in some detail the simplest kind of experiment in which coherence effects 
are important — i.e. Young's (or Michelson's) observation of interference fringes — we can 



^ It is easy to show that 



F{(o) = (l/7c) Re I y"(x) exp[/ cox] dx 
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identify the significance of a light source's degree of first-order temporal coherence and 
demonstrate how that coherence can be measured. Consider the following basic 
experimental configuration: 



Source 




Colli mati ng 

Lens 

Screen 1 

The time integrated/averaged intensity at a given point P is proportional to 




l{P)=2(]^E^^\p,t)]y{E{Pj)E*{P,t)) 



I 



Observation 
Point 

Q 



Screen 2 



[VIA- 17] 



The real instantaneous field the point Q is derivable from the following analytic signal 

E{Q,t) = Q E{P„t- t,) + Q E{P^,t- t^) 

= C^E[P^,t- s^lc) + C^E{P^,t- sjc) ' 



[VIA-18] 



It follows that the time integrated/averaged intensity at the point Q is proportional to 
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l{Q)= QC^(E{p„t- 1,) ^{P„t- {E{p„t- 1) ^{P„t- h)) 

+C,C*{E{P„t- t,) E\P„t- t)) + C,Cl{E{P,j- t,) E*{P„t- t,)) ' 

If the field is stationary , we may with confidence shift the time origin in the various 
expressions so that 



[ VIA- 19 ] 



{E{P„,t-QE*{P„,t-t„))={E{P„,t)E*{P„,t)) 

= /(^„)=ii?(o) 

{E{P„,t- Ef{P^,t- tj)= {E{P„,t+t„-t„) E*{P^,i)) 



nm \ m nf r 



The normalized coherence function 



yS(x)= 



4"(x) 



[ VIA-20a ] 



[ VIA-20b ] 



[ VIA-21 ] 



is often called the complex mutual degree of coherence . Finally, we write the 
general interference law for stationary optical fields as 



i{Q)= liQ) + hiQ) + 2 JTJej y[iM yn 

where I,{Q) = \C,f I (Pj and I,{Q) = \C,f I {P,) . 



^2 ^2 ^ 



\ C J 



[VIA-21a] 
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Thus the degree of first-order temporal coherence may be expressed directly in terms of 
well-defined experimentally observable quantities — viz. 



where I[Q) and I„{Q) respectively, proportional to the intensity measured at Q with 
both pinholes open or with only the nth pinhole open. The visibility of the fringes 



MODELS OF Chaotic Light Sources 

A Collision-Broadened Source - An example of homogeneous- 
broadening: 

Perhaps, the simplest model of a chaotic source is that of a collision interrupted 
oscillator. In this model one pictures an ensemble of atomic radiators each of which 
emits a field of constant amplitude, oscillating at a fixed frequency (co) . The phase of 

the continuous output of each radiator is randomly modulated by the dephasing effect of 
collisions which occur random at some mean collision rate of x~'. The figures on the 

next page shows samples of the field radiated per atom (left-hand graphs) and the 
intensity radiated per atom (right-hand graphs) for ensembles of 50, 100, 200 and 300 
collision interrupt atoms.^ . This model is treated "experimentally' in the discussion 
entitled A "Toy" Model of a Randomly Fluctuating Optical Field which is located at 




[ VIA-21b ] 



given by 




[ VIA-22 ] 



7 



The figures show samples of the sums 



1/ r 

I 2^ sm(|) ,. I / and I 2^ sm(|) ,. I / for atoms with the ph; 





lase 



radiated by each atom distributed uniformly between and 271 . 
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http://www.deas.harvard.edu/courses/ap216/lectures/ls_3/ls3_u6A/toy/toy.html. From 
that discussion, it is reasonable to write for a collision-broadened radiator 

Y®(x)= exp[-/ (co)x- (x/Xo)] [ VIA-23 ] 

(see a graph of this function below). From Equation [ VIA- 15 ], we see that the 
spectrum of a collision-broadened source is given by the Lorentzian form 

itw) = ^l'^- + (V^o)'j [ VIA-24] 

Field and Intensity Fluctuations for Collision- or Doppler-Broadened Sources 



Mean =0.00494765 

Variance/ Mean =0.95 




10 20 30 40 50 
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Doppler-Broadened Sources -- An example of inhomogeneous-Broadening: 

A rather more intricate model of a chaotic source is that of a Doppler-broadening 
system of oscillators. In this model one pictures an ensemble of randomly moving 
atomic radiators each of which emits a field of constant amplitude, oscillating at a fixed 
frequency, but the observed frequency of a given atom is Doppler- shifted with respect 
to that fixed frequency. In the Doppler model, we write the sum of fields emitted by a 
ensemble of moving radiators as 

E{t) = £o Eexp {-/ [o), t-<^ ,]} [ VIA-25 ] 

where the o), 's are the frequencies of individual atoms which are Doppler shifted from 
the mean radiated frequency (co) . In this model, the effects of collisions are neglected 
and, consequently, the first-order correlation function is written as^ 



Since the dynamics of the individual afoms are uncorrelate, all cross-terms vanish. 
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/ 

[VIA-26] 

= E^Y.(exp[-i(ii. x]) 



where x = t^ - Z2lc-t^ + zJ c . The sum may be converted to a integral over a Gaussian 
distribution of Doppler-shifted frequencies and, hence. 



[E\z^t) E{z2t2))= El J exp[-i(o, x] exp|-(co, - ((o))725'j ^/co, 





[ VIA-27 ] 



The corresponding degree of first-order temporal coherence is (see a graph of this 
function below) 

Y® (x) = exp[-/ (to) x-6 ' xV2] [ VIA-28 ] 

and the associated normalized power spectral density follows from Equation [ VIA- 15 ] 
as 

i^co) = [inh']"' exi[-{co- <co))725'j [ viA-29 ] 
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Degree of First-Order Coherence - Absolute Value 




Spatial Coherence: The van Cittert-Zernike 

In the rudimentary treatment of the Young interference experiment recapitulated above, we 
assume an ideal point source and, thus, neglect any effects of spatial coherence . In 
particular, we assume that the field at points Pj and P2 due to a given radiator are in time 

synchronism. However, when we deal with extended sources, we must enlarge our notion 
of coherence The subject of spatial coherence was first developed by van Cittert and, later, 
it was extended by Zernike.^ To define the problem, let us consider the following 
geometry: 



9 See P. H. van Cittert, Physica, 1, 201 (1934) and F. Zernike, Physica, 5, 785 (1938). 
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The task before us is to deteraiine the so called mutual intensity /(P^Pj) due to the 
extended source which might be observed for the two points located on the observation 
plane. For simplicity, we assume that extended source is in a plane parallel to the 
observation plane and that it can be divided into cells of statistically independent radiators. 
Thus, the total field at any point on the observation plane is given by summing the fields 
due to each of the cells ~ viz. 

E{P,t) = Y,E,{P,t) . [VIA-30] 
The mutual intensity is then defined as 

= l.'\UPvt) El{P2,t)) + Il,{E,{P„t)El.{P„t)) ■ [^^-31^] 

k k'*k k 
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Because of the assumed statistical independence of the cells, the cross-terms vanish so that 

■/(^,^2)=E(4(^P^) KiPi^t)) . [VIA-31b] 



It the spirit of the Huygens principle , we assert that component fields radiated by each 
coherent cell are spherical waves of the form 



E,{P,t) = E 



PS,. 



[expj^-/(a)) t-'-^ 



V c jj 



PS, 



\. [VIA-32] 



where E[Si^,t) is field radiated by each element of the source and P 5^ is, obviously, the 
distance between the source and observation positions. Therefore, we can write 



[E,{P„t)E:{P„t))=(^E 



c )/ [ [P,S,][P,S,] j 



/ / ^ */ , ^\\^Mi{<^){PiS,-P2S,)/c] 



[ VIA-33 ] 



In most instances, the path difference will be small compared to the coherence length of the 
cell, so that, in the limit of many small cells, we obtain the famous van Cittert-Zernike 
theorem in the form 



7(5)- 



Sou 



[n s] [p^ s] 



dA^ [VIA-34] 



and the first-order degree of mutual (spatial) coherence is defined as 
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, , J(P,,P,) 



" * Sniirre 



where 



Source 



2R 

and, hence, we may approximate Equation [ VIA-35a ] as 



[ VIA-35a ] 



nP„)= II r^^^. • [VIA-35b] 



For most problems of interest, we may take 

P. S = ^R'Hx„ - -^y-^R+ ^'"- " [VIA-36] 



- . Source 

Yn — 



[ L Jj 

||/(^,Tl)rf^jTl 



[ VIA-37a ] 



Sou 



Further, we see that in the Fraunhofer or far-field approximation, the mutual 
coherence is the Fourier transform of the source intensity distribution — viz- 
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I {^,^)exp[i{k)^{x,-x,)/R] exp[ i{k) ^{y- y,)j R] d^ 



^ Wee ! [ VIA-37b ] 



Source 

— and, consequently, the measurement of the degree of mutual coherence of a 
remote source provides a means to measure the angular size of that 
extended sources -- i.e. stellar objects. 

We see then that a complete first-order description of the coherence properties of a classical 
field is embodied in the complex first-order degree of spatial-temporal 
coherence which is defined as 



f'^ {\t,,\t,)^ 1^ " [VIA-38a] 



where 



'\E*(\,t)E(\,t,)) 
(|£(r„rjf\/|£(r„r,)f 



1 

{e\\, t) E ^ j E%,t, + t) E(i,,t,^t)dt ■ [ VIA-38b ] 



Free Space Propagation of Coherence Functions: 

We derive here an equation which allows us to generalizes the notions imphcit in the 
van Cittert-Zernike theorem by providing us a description of how coherence propagates 
through space. We start by assuming that the analytic signal representing the field of 
interest satisfies an wave equation of the form 
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£ (?„ , = A ^ £ (I ,t„) [ VIA-39 ] 

c at.. 



If we multiply through by E (f,, fj) , we see that 

V',r{i,t„i,t,)=V',{E(i„t)E%,t,)) 



= -J^^,E(i,,t,)E%,t,)y [VIA-40a] 



c 
1 



c dt^ 
Similarly 

vl r (?, t„i,t,) = -^,^,r (?, ?2 ^2) . [ viA-40b ] 

Therefore the coherence function must satisfy the fourth-order equation 

1 

Vi v^r(f„f2;x)= -,—,r(i„i,;x) [via-41] 



The imphcation of Equations [ VIA-40 ] and [ VIA-41 ] is that coherence is a field that 
propagates through space and may be treat by methods which apply to other field 
variables. In fact, the content of the van Cittert-Zernike theorem is a description of how the 
coherence of extended source propagates to remote observation points. This point of view 
is particularly valuable in the treatment of wave propagation through random media. As a 
beam propagates through such a medium, the effect of random scattering events scrambles 
the phase of the beam and the simple notion of a field rapidly loses meaning. However, the 
evolution of the coherence function provides an important measure of the statistical 
properties of the random medium. 
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Higher-Order Correlation Functions -- A classical explication of 

THE FAMOUS HANBURY BROWN-TWISS EXPERIMENT: 

Generalizing Equation [ VIA-38a ], the degree of nth-order spatial-temporal 
coherence can be defined as'" 



j(i£(f.. 4) (f- 'jf ) (1^ (f- '-.)f) (f.-'=.)r) 



[ VIA-42 ] 



In particular, the degree of second-order temporal coherence is defined as 

/ \F.(r , t\\^\F(r , t+T\^ \ I 



(2) — X ■ I \ I 

^ ^'■-f'^^^ 77^ = 7^7- u2 [VIA-43] 

We shall see that this function is an important measure of the relative timing of of intensity 
fluctuations.' ' As a first step, let first find an expression for the average intensity radiated 
by a collection of independent oscillators — viz. 



/|£(r,,,r)| |£(r,,„f+x)| \ (/(r,e,,r)7(4f,^+x)') 



l'' R. J. Glauber, in Quantum Optics and Electronics, Les Houches, 1964 (edited by C. DeWitt, A. Blandin, and 
C. Cohen-Tannoudji), p63, Gordon and Breach (1965). 

11 By the famous Schwartz inequality 



<l\.\'\/\uA 



and quite generally 



so that 



Y*" (f.ef.^)^l 
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I - (/ {t))= ( |^(0I> Is^E^t) E;{t)) 

\u I [VIA-44] 

=[Boc/2]i:(|£,.(r)f) 

i 

In this expression all of the cross-terms are presumed to vanish, since the radiation from 
any given atom is statistically uncorrected with that of any other atom. For the coUision- 
and Doppler-broadening models discussed above, the oscillators differ only in phase and 
this average can be expressed as 

I = (/ [t]) = A^,,„ [8„ c/2] \af. [ VIA-45] 

To obtain the root-mean-square deviation in the cycle average of the intensity, we first 
calculate 

([/ {')])= [Kcnr { \E{,)\'y[e,c/2r (zE,m')e:i')e:i') ) 

= [e.c/2l={l(|£,(»)|>I(|£,(')l>,(')r)} 
= [e.c/2]4l(|£,(,)r)+I(|£,(r)r)(|£,(,)r)} 

I i M J 

= hc/2Y |iV,„^( K(^)r) + 2iV^,„„(iV^.„„-l) [( \E,{t)( 
Again, for the collision- and Doppler-broadening models, this average can be expressed as 



[ VIA-46 ] 
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= 2 



[ VIA-47 ] 



and, thus, we find to a very good approximation that 



P {')]')-{{!{>)) 



= I 



[ VIA-48 ] 



which is consistent with the experimental results presented above. Let us now look at the 
two-time average of the intensity from a collection of statistically independent 
oscillators — viz. 



(/(rJ/(r,))=Kc/2]^(|£(rJ|V(r,)f) 

= Kc/2]'/ Y.Eit)E.{t,)E:{t)E:{t,)) 



[ VIA-49 ] 



As a result of statistical independence, most of the cross-terms vanish so that 



(7(r,)/(o)=Kc/2]Ni:(|£,(on£,.(gf) 



+ 



[ VIA-50 ] 



Again, for the collision- and Doppler-broadening models 
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(7(.07(O)=[E„c/2]NA^,,,(|£,.(rf|£,(r,)f) 



+ A^tex.(A^.ext-l) 



[VIA-51] 



and, hence. 



Y'^'(T) = l + |Y«(T)f 



[ VIA-52 ] 



which is important result which holds for any chaotic source — see the figure on next page. 
Notice that this result predicts a bunching of the intensity variations. 




THE H ANBURY BROWN - TWIST EXPERIMENT 
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.Schema'tic diagram of the Hanbury Brovn-Tviss experiment 
demonstrating coherent interference in light from incoherent 
sources. Sources 1 and 2 are separate by distance d1 2^ and emit 
photons of identical wavelengths which are detected by detectors 
A and E located a distance dAE apart. The distance between the 
source plane and the detector plane L. A coincidence betwe&n 
detector outputs results in a composite signal exhibiting an 
-interference effect depending on both d1 2 and dAE. 
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Light 
Source 




n 



Electronic 
Correlator 



The Hanbury Brown - Twist Spectrometer is designed to measurei^ 

([A(fref. - a] [^(fref.^l) " h]} => i{{jfircfA) h{ircf' ^l))" ^ '} [ VIA-53 ] 



This correlation may then be written 

([A(fref>0-A][4(fref>^2)-4]) 



IJ2 



[ VIA-54 ] 



Generally, the results of these experiments are consistent with Equation [ VIA-52 ], but 
some of the more detailed observations require a quantum mechanical interpretation. 



12 R. Hanbury Brown and R. Q. Twiss, Proc. Roy. Soc, A 242, 300 (1957). 
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VII. QUANTUM MECHANICAL LANGEVIN EQUATIONS 38 
A Rudimentary Resevoir Problem :394o 

Consider a complete system which consists of a single simple hannonic oscillator (the 
system component of "interest") coupled to a reservoir of many simple harmonic 
oscillators (the system component treated statistically) 



"The System" 
A single simple 
harmonic oscillator 



^nt 



eraction 




Reservoir 
of 

many harmonic 
oscillators 



In particular, the complete interacting system has a simple Hamiltonian of the form 

^ = ^ laser + ^ res+ ^ int 

= hSla^a+^h(Oj b]bj+^h[gj b. + g* b]a\ 



[VII-1] 



where a is the system variable of interest (the "state of the system") and the b-^ are the 
other system variables (the "reservoir states"). System operators commute with reservoir 
operators at given time, we may easily generate the following set of equations of motion: 



Much of what follows draws heavily on material in the Arizona Books — i.e. 

1. M. Sargent III, M. O. Scully and W. E. Lamb, Jr., Laser Physics, Addison- Wesley (1974) 

2. P. Meystre and M. Sargent III, Elements of Quantum Optics, Springer- Verlag (1992) 

3. Weng W. Chow, Stephan W. Koch and Murray Sargent III, Semiconductor-Laser Physics, 
Springer- Verlag (1994) 

This model is particular valuable in treating electromagnetic intermode coupling problems, but it also provides 
guidence in treating more general problems which incorporate atomic states. 
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a{t) = i [fH,a{t)]/h= -iQ. a{t) - i^^gj bp) [ VII-2a ] 

bj{t) = i [i^,bj{t)]/h= -/CO. bj{t)- igj a{t) [ VII-2b ] 

Fonnally integrating Equation [ VII-2b ], we obtain '*i 

bj{ t) = bj{ to) exp[- / co^. {t - ?(,)] - i g* J dt' a{t') exp[- i(ii.[t- [ VII-3 ] 

The first term in this equation describes the free evolution of the reservoir states in the 
absence of any interaction with the system and the second gives the modification of this 
free evolution as a consequence of the coupling to the system. Inserting Equation [ Vn-3 ] 
into Equation [ Vn-2a ] we obtain 

d (f) = - j 12 a (r) - ^1 g,. f J d t'a [f) exp[-/ (i). [t- 

' [ VII- 

- i 8j bj{t^) exp[-i (0. (r- rj] 

j 

Here the second summation describes a source of fluctuations arising from the free 
evolution of the reservoir states and the first is a feedback through the reservoir which 
might be described as a radiation reaction. To remove the rapid variation in the system 
variable we transform to the Heisenberg interaction picture ~ i.e. we take 

a{t) = A {t) exp[-ii2 d [ VII-5 ] 

As is pointed out in Section 14-3 of P. Meystre and M. Sargent III, Elements of Quantum Optics, Springer- 
Verlag (1992), the statement that "System operators commute with reservoir operators at given time. .." is a 
very critical point. From Equation [ 11-3 ], we see that as time evolves, the reservoir operators acquire some of 
the characteristics of the system operator. Thus, it is crucial in this development to maintain a strict and 
consistent ordering of operators at different times! 

The transformation does not alter the system's communtation relationship, so, of necessity, 
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so that A [t] = g. ^ \ ' dt A{t') exp[- j (o). -O) {t-t)\ + F{t) [ VII-6 ] 

where ^(0= 8j ^,(0 exp[i(i2-(a .) {t- t,)] [ VII-7 ] 

i 

is the so called noise operator. 

The first term on the RHS of Equation [ VII-6 ] may be reordered as 

Y}gj r J ' dt' A{t')^xA-i[(Si-Q) {t- t')] 

i 

= \' dfA{f) XI 8j r exF[-/ (oij -a) {t- t')] 



[ VII-8a ] 



By assumption, the frequencies of the reservoir oscillators are closely spaced and widely 

distribution so that the function X| gj | exp|-/ (to^. [t- is very sharply peaked at 

j 

t'= t and has a width of the order of the inverse of the reservoir's frequency bandwidth. 
Thus, A[t') may taken out of the integration and the limits of the integration may be 

extended to infinite - i.e. 



Yj[gj r [ rf?'A(0exp[-/^-t2) [t- t)] 

i 

W LI Sj f \^dt'ex^-i ((Hj-Cl) {t- f)] 



[ VII-8b ] 



Again following arguments explicated by Heitler,'^^ jji^t 



In Chapter 11, Section 8 of W. Heitler, The Quantum Theory of Radiation (3rd edition) 
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J dt exp[-i(i2-co .) [t- t)\ =>7t 5(^2-03 .)- Tr —!— [ VII-9 ] 

~°° ~~ CO j 

If we neglect the imaginary part of Equation [ VII-9 ], which leads to a Lamb-t5^e 
frequency shift. Equation [ VII-6 ] leads directly to a quantum mechanical version of the 
classical Langevin equation ^4 — viz. 

A{t) = -^A{t) + F{t) [VII-10] 

where Y= 2% p{Q)\g{a)\\ [ VII-11 ] 

is the quantum mechanical Langevin drift coefficient. In optical applications the average 
intensity or number operator — i.e. (a'(?)A(?)^ — is the experimental quantity of 

greatest interest. A useful equation of motion for that operator may obtained by making use 
the Langevin equation ~ i.e. Equation [ VII-IO ]. To that end we first write 

|(aV) A (r))= -y(aV) A(r)) + (FV) A (0)+(aV) F[t)) [ VII-12 ] 
To obtain {F\t) A [if) and [A\t) F {t)) we make use of the identity 



^ In the classical Langevin theory of Browian motion, the equation of motion for drift velocity of a particle 
suspended in a liquid is given by 

mu (r) =-m Yu [t) + Fj^t) 
where F is a damping constant and F^^i^t) is a random, rapidly fluctuation force with zero mean value. In this 
classical problem the famous fluctuation-dissipation theorem is expressed as 

T=^D^^ = ^[^dx{FXO)FXx)) 
where D^^ is the so called Fokker-Planck diffusion coefficient. 
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A(t)=A{t-At)+\' dt'A(t) [VII-13] 

so that 

'F\t) A (t)) = (F\t) A (t- At))+ \ ' dt' lF\i) A {A) 

[VII-14] 

= {l^{t)A{t- A.))- 1 ]]jt'{P{t) A(0)+ ]'jt'{F^{t)F{t)) 

The first term on the RHS vanishes since cause cannot precede effect. That is, M^i) 
cannot be correlated with a future value of the fluctuating force. Similarly, in the second 
term {F\i) A[t'^f vanishes except in the infinitesimal small inteval where t' = t. 

Therefore 

[F\ t)A{t))=\' Jt' {f\ t)F{t)). [ VII- 15a ] 

If the random force is stationary in time 

[F\t) A{t))=Y^ dx{F\t) F{t+x)) [ VII-15b ] 

and, finally, if A? is long compared to the random force correlation time 

{F\t) A {t)) = I" dx {F\t) F {t+x)) = ^ J^°°Jx {f\0) F (x)) . [ VII-15c ] 

Therefore Equation [ VII- 12 ], the equation of motion for the average number operator, 
becomes 

I {A\t)A{t))=-y {AXt)A{t))+ J Jx {f\0),F{x)) [ VII-16 ] 
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which is one form of the famous fluctuation-dissipation theorem. 

If the reservoir is in thermal equilibrium, we may directly evaluate the noise operator 
correlation to wit. 

(^V) ^(0),,™= E 8* & (^;(0) b,{0)l^^^ exp[-/(t2-co,) t] exp[/(n-a),) /] [ VII- 

Since the operators of the reservoir commute we can write 
(bmbJO)) =(n) 5,., 

\ J\ / t\ / 1 ffigy^ \ J I therm 

= ^ In^ Lexp^-n . p ^coj \ [ VII-18 ] 

= {exp[p;za)J-l}"5,., 

and Equation [ VII-l? ] becomes 

[H') F (<-)),„.= I f ( n. exp[-,-(n-<0,) (,- 0] 

i 

= 27r|^a)fp(t2) (n(t2)X^_5(r-0 

= ^'\<^)),,eJ^'-') [VII-19] 

\ I therm ^ ' 

/ \ y 

where ( D , ) = — ( n (C^ ) is of the form of a diffusion coefficient in the 

\ / Iherm 4 therm 

classical Langevin theory. Two-time correlation functions of this form are characteristic of 
Markoffian random processes and represent fluctuations in a reservoir with essentially 
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zero memory. Substitution of this Markoffian correlation function into Equation [ Vn-16 ] 
yields an Einstein relation of the form 



General Resevoir Problem 

As the previous section demonstrates, the effects of random fluctuations in an assemblage 
of harmonic oscillators may be accounted for in precise detail. The problem can be 
analyzed completely since the first term on the RHS of Equation [ Vn-6 ] ~ the damping 
term — contains only the system variable. In general, this feedback term will included 
reservoir variables as well and the simplified model breaks down. The general case may, 
however, be analyzed if we take the simpler case as a guide. In particular, let us suppose 
that a set of systems operators { } are coupled to the reservoir and satisfy a set 

Langevin equations of motion as 



\ I therm df^ 



I thi 



'.erm 



[ VII-20 ] 



[ VII-21 ] 



where D^[t) is the drift term and i<J (r) is the Markoffian noise operator appropriate to the 
system operator A^[t)- From the earlier results, we are, thus, assuming that the two-time 

correlation of the random fluctuating force is of the form 



{^^(0^v(0)=2^.v 5(^-0 [VII-22] 



Again using the identity 



Mt)=Mt-At)+\ljt'\{t') 



[ VII-23 ] 



we find 
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'M)Fit))={\{t-At)FXt)) 

r, r, [VII-24] 

As discussed in connection with Equation [ VII- 14 ], the first and second term on the RHS 
in Equation [ VII-24 ] vanishes since cause cannot precede effect. Therefore 

(A^{t)Fit))=\'jt(F^{t)Fit)) [VII-25] 

and by Equation [ VII-22 ] 

{Mt)FXt)}=D^^ [VII-26] 
We may then form the following equation of motion 

= {D^ A\+{F^ A\ + {A^ D\ + {\ F\ 
which leads to a generalized Einstein relationship 

2^,v = ^(A^ A,)-{D^ D,) [VII-28] 

With this equation it is possible to calculate the diffusion coefficients from the drift 
coefficients, provided one can independently calculate the equation of motion for (^A^^ A^y 

We can also use Equation [ Vn-21 ] to calculate the expectation values of two-time 
correlation functions. In paryicular, if multiply it by A^(?') where and average we 

find 
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jjMt) ^v(O) = {DM Ait)y[F^ [t] AXt)) [ viI-29 ] 

Since in the Markoff approximation, the present system operator Aj^^t') cannot know about 
the future noise source ij^ (?) , (0 ^v(O) vanishes and 

^v(O) = {DM ^v(0) [ VII-30 ] 

which is the so called quantum regression theorm which states that the two-time correlation 
function \\{t) ^v(O) satisfies the same equation of motion as the single-time (?)) • 
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Semiconductor Langevin Equations 

As discussed earlier, the effective Hamiltonian for optical interactions in a semiconductor 
may be written as 



^ eff - 



Thus, the equation of motion motion for the effective dipole operator ^ is given by 



d i r -I 

-k k fcL eff -k kj 



=-(y+ i CO-J 6_^ci-^+ig-^a(cil + 6 t B 1) + /, ^„ 



[ VII-30 ] 



where y is the damping, a)- = £ H —-\ —, and fs -^- is the k-dependent noise 

2m 2m. ' 

ec hv 

operator. We transform the various operators to a frame rotating at the electromagnetic 
field frequency (0^ — viz. 



a{t) = A{t) exp[-a), t] 

fe-va,= U=^^, exp[-co, t] 



[ VII-31a] 
[ VII-31b] 
[ VII-Blc] 



and find the transformed equation of motion in the rotating frame ~ i.e. the Heisenberg 
interaction picture ~ as 

^ a- =- ( Y+ / CO- - / CO,) a- + /g - A (a 1 a- + ^ t ^ _- - 1) + F^^ [ VII-32 ] 
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Applying the generalized Einstein relation ~ i.e. Equation [ ¥11-28 ] ~ to the operator 



we obtain 



Similarly, 



33 



[ VII-33a ] 



[ VII-33b ] 



where al,= (l-al a-^(l - b l^b _-) . 

The equation of motion for the laser field operator in a frame rotating at the laser frequency 
co^ may be written as 



Kt) - iHg^ + [ VII-34 ] 



where Q. is the cold cavity resonant frequency and where we have replace the damping 
parameter y in equation [ VII- 10 ] by the experimental parameter (ajQ . If we consider the 
Equation [ VII-32 ] in the quasi-equilbrium limit, we find 



O-lt) = [y+ i (D,- /CO,]"' [ig^ A{t)(ala^ + BXB _^-l) + F^^{t) 
= !Z)(a)- -CO,; y) f /g- A{t) [0% - a^) + 



[ VII-35 ] 
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where is the complex Lorentzian demomnator.'*^ Thus in the quasi- 

equilbrium Umit, Equation [ VII-34 ] becomes 



r 



1 



^^^""[^^'^"""'^^"^I^^|'^"^""^'^H^'3-<7oo)j^(0 + ^a(0 [VII-36] 



where 



FA{t)= F{t)- 2)(co--co,; y) F^_^{t) 



[ VII-37 ] 



Following arguments similar to those used in deriving Equation [ Vn-16 ], we easily derive 
an equation of motion for the laser photon number operator ~ viz. 



^^-l Lkf 4^i-(xi/,y){c,,-cj'^{A\t)A{t)) 
^ ^ i J [VII-38a] 



+ \^dx{F}{0),F^{x)) 



or 



-{A\t)A{t))=-\^-{R,,-R^) 



loo 

la 



(A'(t) A(t))+ {fI{0).F,{t:)) [ VII-38b ] 



where 



' k 



[ VII-39a ] 



It may be recalled that in the notes entitled Lasers: Models and Theories for convenience we indroduced the 
Lorentzian functions - viz. 

2^ M — v; w) = [ I ( M — v) + w] the complex Lorentzian denominator 

v; w) = |^(m- v) the dimensionless Lorentzian function. 
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and ^00 = I Lkf r; y) (5 00 . [ VII-39b ] 



Y 



k 



Assuming that the cavity damping and carrier scattering reservoirs are uncorrelated, the 
two-time correlation function may be written as 

k k' 

and if the noise operators for different k are uncorrelated 

{Plit) m) = {F\t) Ikl x(co,-a).; y) (i^^(r) . [ VII-40b ] 

' is 

Finally, using Equations [ VII- 19 ] and [ Vn-33a ] we obtain 

(^V)^a(0/= I W«)/,,.™+^Eld'4«i.-«.;Y)cy33js(^^ . [VII-40C] 

and, thus, the " A ' A " diffusion coefficient for the laser field is 

CO / / x\ 

2^.u = -(n(co))^^_+i?33. [VII-41a] 
By a similar argument, the " AA ' " diffusion coefficient is 

^^AA^ = I [("N)...™+l] + ^oo . [ VII-41b ] 



Using Equation [ VII-41a ], Equation [ VII-38b ] becomes 
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CO 



g ~ (^33 ~ -^Oo) 



(A*(»)A(,))+^ («(«.)>,.„.+ «: 



33 



[ VII-42a ] 



Similarly 



(A(.)AV))+|[(n(a))) 



therm J 



[ VII-42b ] 



As check on the consistency of these result, we see that the form of Equations [ VII-42a ] 
and [ VII-42b ] guarantees that the commutator \A{t), A^(?)] is independent of time ~ viz. 



dt 



'iKt)At)\)= 



(0 



— - (i?33 - R^) 



I^A{t),A\t)\)+^ + {R,,-R,,)= [ VII-43 ] 



to incorporate saturation effects into the noise spectrum we need to write a Langevin 
equation for the carrier-density operator ~ viz. 



— *k *k~~['^eff'^k^k]''" {damping, pumping, and noise terms} 



dt 



[ viI-44 ] 



where A*'^ [l ~ ^ ^ t] is the pumping due injected carriers, al a-^ B }^ B is the radiative 
recombination (spontaneous emission) term, Y„r ^ « ^ is the nonradiative recombination 



term, 



is the relaxation due carrier-carrier scattering, and is the 
carrier noise operator. In the quasi-equilibrium approximation ~ i.e. from Equation 
[ VII-35 ] — we may sum Equation III -44 ] to obtain a Langevin equation for the total 
carrier-density operator 
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A AT = A- 1 y r aU-B'-B - - Y„, 

y L-i k k k -k -k 'nr 

lA'A 2 - - [VII-45a] 

- — TT &kf 4^k-«.; y) ((^33 - 

' k 

where 

^.(0= + 40 2^co--co,;y) A(r) lD>t-co,;Y) [ VII-45b] 

k 

NOISE SPECTRA: 

To calculate noise spectra we make use of semiclassical transformations ~ viz. 

= (£(0/^»J e'^rf-'^lO] [ VII-46a ] 

A'(r) = (4r)/2;»J exp[j>(r)] [ vil-46b ] 

£(?) = 'E„^ [VII-46C] 

(t<?) = ln[A(?)/AV)] [ VII-46d ] 
where T.^^ = [h (ajs-o V] is the electric field per photon. Using the expression 

=yD [VII-47] 

we see that 
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a_E a£ dE dE 



dA dA^ 



4E' 



AA' D^,^+A'A D^^] ^^[d^,^ + D^^] 



1 1 r 1 



[ VII-48a ] 



[ VII-48a ] 
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VIII. NONLINEAR OPTICS -- QUANTUM PlCTURE:^^ 

A Quantum Mechanical View of the Basics of Nonlinear Optics 

In what follows we draw on the discussion of the density operator in Review of Basic 
Quantum Mechanics: Dynamic Behavior of Quantum Systems, Section II of the lecture 
set entitled The Interaction of Radiation and Matter: Semiclassical Theory 
(hereafter referred to as IRM:ST). The macroscopic polarization is given by 

(p)=Tr{pP) [VIII-1] 

We take the total Hamiltonian of a particular system in the form 

!H = i^o + ^,„.+^™,do. [VIII-2] 

where f)-f random is ^ Hamiltonian describing the random perturbations on the system by 
the thermal reservoir surrounding the system. Thus 



where 

dt 



[ VIII-3 ] 

i"ela;( 



= ^[p,^.a„doJ [VIII-4] 

relax 



To find the nonlinear susceptibility, we make use of the following perturbation 
expansions: 

p= p("'+pW+p<2) + ... [VIII-5a] 



See Nonlinear Optics — Classical Picture which is Section VII in the lecture set entitled On Classical 
Electromagnetic Fields (OCEF). 

See, for example. Chapter 2 in Y. R. Shen's Principles of Nonlinear Optics, Wiley (1984). 
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with 



/p\_/p(0) \ ,/p(l) \ ,/p(2) \ , , 



/p{a) \=Tr(p*"' P) 



[ vni-5b ] 
[ vni-5c ] 



By substituting these expansions into Equation [ VIII-3 ] and equating terms of like 
order in ^ we obtain the following hierarchy of equations: 



(0) 



x(0) 



(1) 



ap^ 



dt h 

-;{[p»,^,].[p'»u.]}.f 



relax 

ap*' 



ap 



(3) 



dt 



relax 



relax 



relax 



[ Vni-6a ] 



[ vni-6b ] 



[ vni-6c ] 



[ VIII-6d ] 



Following earlier considerations, it is reasonable to write 



l«>-r„„-(Hp|«>-r„„'P„ 



relax 



[ VIII-7 ] 



Since the perturbing field is resolvable into an appropriate set of Fourier components 
(either discrete or continuous), we may write 
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^int = intK) exp(-/a),?) [ Vni-8a ] 

k 

P = Lp ((0 J exp (- / [ Vni-8b ] 

k 

and resolve Equations [ VIII-6 ] into a hierarchy of algebraic equations. The first 
member of that hierarchy becomes 

9l\^k)[i{^nn-^k)+r.n]={{n\{9'''^ 

^ [ Vni-9a ] 



or^« 



p!i-K) = ^(H^.mK)k){pl?-pit }2^K„-co,;r„„.) . [viii-9b] 



The second member of the hierarchy becomes 

PwK){ ^■K'-«.)+r„„,} = ^X<'^|[p»(a),,),i^.,(a),-a),.) \\n') [Vni-lOa] 



or 



p2'K) = ^ 2>K„^-co,;r„„,)H{ pli„(to,,) <n"| ,,(a),-Q),,) 

^ « " [Vlll-lOb] 

-(n|5/,,K-(0,.)|nOP?vK)} 



we presume in this development that ^ is a Hermitian operator. 

Recall that M — v; w) = [t ( M — v) + w] is the so called complex Lorentzian denominator. 
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Thus, we see that this expansion yields results akin to those embodied in Equations 
[Vn-8 ] through [ VIII- 13 ] in OCEF. 

Two-PHoton Absorption 

Let us consider two-photon absorption by a single atom (or by a set of paired k -states 
in a semiconductor as discussed Section 7 of IRM:ST). 



COb 


\C02 

— ^ 







^\co^ 




COb ^ 


COa 


^/^COi 


^ 



Two-Photon Absorption Process 

Tfrom earher discussions, the appropriate interaction Hamiltonian is given by 

H,, =-e f - E<^' (r,?) =-e ? [e^' (R,r) + E<^' (r,?)] [ VIII-U ] 

In the appendix to this section we establish — viz., in Equation [ VIIIA-12 ] — that the 
second approximation to the transition rate is given by 



co,-co,„ 



6(ol),. -cd^) 



Thus for the two-photon process 
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J?) 



2%e' 



</| ? |m)(m|r E<"'|/) 



5((0,.-(0^) [VIII- 12] 



For two beam, two-photon processes 



2 Tie' 



IE 



(/|r.Er'|m)<m|r.Er)|/) 



(0,-0), 



+ 



(/ir-ErvxHf-Er'io 



[ VIII- 13a ] 



5(0)1 +0)2 -(Oy^) 



and for one beam, two-photon processes 



2 Tie' 



r(2) 



E 

/ 



</|rE(^'|m)(m|rE<^'|/) 



0), -(0, 



S{2a),_-a),)[VIII-13b] 



If we use once again the factorization used in Section VI of this lecture set — viz- we 
write the initial state as | ? ) = |i = {''''"})~| (0 )| ^i^'l ) where | A(/') ^ and 
I F[i"] ^ are, respectively, the initial electronic (atomic) and field (photon) states. 

Using closure on the intermediate photon states, we find for two beam, two- 
photon processes 



xTr [pE'-^E^-' E^^'E*^' ] 



[ VIII- 14a ] 



where 
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^ (1,2) ^£ 



■(A(/)|r|A (m)) e,e, (A(m)|r|A(/)) 



+ 



(A(/)|?|A(m))e,e,-(A(m)|?|A(/)) 



CO. -co 

2 m 



[VIII- 15a] 



and for one beam, two-photon processes 

1 2 Tie' 



(1,1)1 5(2co-co,) 

Tr [ p Ej-' E|-' E{"' eW ] 



[VIII- 14b] 



X 



where 



M (1,1) 

m 

Four Wave Mixing: 



<A(/)|r|A(m)) e, e, ■ (A(m)| r |A (/) ) 



COi-CO„, 



[VIII- 15b] 



Four Wave Mixing Spectroscopy: 

We explore here a particular set of four wave interactions. To that end, we consider an 
input field which consists of two plane waves 



E«(r,f)= E„(ki,C0i]exp[/(ki r-C0i t)] 

+E„(k2' C02)exp[/{k2-?-C0 2f) 



[ vni-16] 



-i-c.c. 



where frequencies cOj and cOj lie in the visible part of the spectrum. As Equation 
[ VIIT3 ] of OCEF informs us, this input incident on a third-order nonlinear material, 
will directly generate a component of nonlinear polarization at a frequency 2C0, - CO 2 

- viz. 
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(3)1 



X 



exp (2k[- kj) ■ r| exp [-/(2CO1-CO2) 



[ Vni-17] 



which may be associated with the destruction of two photons at (Oj and the creation of 
photons at CO2 and 2a)j-C0 2. 

L 




Four wave mixing - direct process 

If, however, the material lacks a center of symmetry there is also an alternate path to 
obtain 203,-0)2 radiation. In a second-order material [ VIII-3 ] of OCEF informs us 
we can also have a nonlinear polarization at 200,-0) 2 



<^'pnf,r)={2xS;(o),-o)2;a),)Ep(k,-k2. 03,-0)2) E,(k„o), 

+xS;(203,;-0)2) Ep(2k, 2o),) E;(k2, 0),)} [ Vni-18 ] 

xexp[i (2ki- £2) r] exp [-/ (203,-0)2) ^] 
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if we have driving fields Ep[k,- k,. cOi-WjJ and E|3[2kj, 2a)iJ. Such fields are, in 
fact, generate by the input field through the nonlinear polarization components 

<^'prtk,-k2.(0,-(02)=2%S;((0,;-a)2)Ep(k,a)jE;(k,,a),) [VIII-19a] 



and 

<^'pr{2 20),) = 2xil{(0-(0,] Ep(k„a),) E^(k, o)J [ VIII-19b ] 




Four wave mixing - indirect processes 



We may solve for Ep[k, - kj.co, -COj] and Ep[2k,. 2(oJ from Equation [ VIII-23 ] in OCEF — viz- 

{|Qf Sap - QaQp - {(oVc'eJ 8„p (o))} Ep (q, (o) (o' Pa'^^'lO. co) [ VIII-20a ] 
or usmg the operator 3(k) = (k ■ k) 1 - k k = |kf |l- k k} defined in OCEF 

{ 8-„; (co)3ps (q) - {co7c>„s } (q, (o) = ((o'/c') e^; (co) Pp<'^"^{Q, (o) [ VIII-20b ] 
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If we define 

g„s{Q;a))^{8o8-„;(a))3p4Q)-{a)7c')6„8 } [Vni-21] 

we can, in principle , write 

E Jq, 0))= [coVc^] g (Q ; (o) z^l (co) P8*'''{Q. O)) [ VIII-22 ] 

Therefore, for a material lacking a center of symmetry, the complete nonlinear 
polarization may be written 

V:\f,t)= 3exp[/(2k - k,)-f] exp[-/(2co -0),) t] Ep(k,co) E^(k,(oJ Ejk,. co^) 

x{xaM K;wi;-«2) 

4 [VIII-23] 



3 

2 



+ 



Notice the "resonance" in the second term when |k,- kj] and COi-COj are "tuned" to the wave number and frequency 
of coupled electromagnetic excitations or normal modes such as plasmons, polaritons and magnons. 



Degenerate Four Wave Mixing (DFWM) - Phase Conjugation: 

Consider now the reflected waves generated by the direct process — i.e.. Equation 
[ VIII- 17 ] — when the input field is the sum of two counter propagating "pump" fields, 
with [kp,(o] and [-k^^,co] , plus a complex wave E,„(r, i) represented a range of wave 

numbers components ||kf},coj — viz. 
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E„(f , t) = |e1:' exp [/ (k^, • f -(^,)] + E<-' exp [/ {-k,, • f-^ _) 
+^e|^'* exp[z(k^ -f-^l^/)] ^exp[-zO) t] + c.c 



as shown below 



input 
signal 





/ 

/ nonlinear 
' crystal 



[ Vni-24 ] 



According to Equation [ VIIT17 ], reflected waves are generated by nonlinear 



polarizations components with |^|kf},3tt)j, |{2 k^ + kf},OL)j , ||2 k^- kf},-co 



, and 



. The latter component is the so-called phase conjugate term and may 



[{kj ,-(0 
be written 

*^^Pr(f,0= 3%L;\5(a);a);a))E<-' E«(^, a),)exp[-/((^,+(^_)] 



J^Eg'' exp[-z(k^ r-(|)f)] Uxp[-/co fj + c.c. 



[ VIII-25 ] 



which may be interpreted as the time reversed image of E,„(r, t) . 



TIME REVERSED IMAGING 
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i ncident signal 




The Raman Effect 

The Raman Nonlinearity: 

In our discussion of the dielectric susceptibility thus far we have been implicitly 
assuming that the nuclei of system were clamped into fixed positions and have focused 
on the nonlinearities associate with electron dynamics. This is too limited a view. In 
general, we may expect that polarizability of, say, a molecule is a function of the 
nuclear positions and that the dipole moment can expressed 



|u I denotes a collection of nuclear positions where u is the displacement of the s 

nucleus from its nominal equilibrium position. Further, we may expand these terms in 
powers of the displacement — viz. 



P„({fi'1.E)=pr({s"'}) + <..p({0"'}) 




[ vni-26 ] 





(0) 



1 




a 



[ VIII-27a ] 
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• S({a'''},E).a„J{a«})Ep==a<;'Ep+t4iEpui'' + .-^ [ VIII- 27b ] 

s-l 



The second term in Equation [ VIII-27b ] is the so called Raman nonlinearity, 

has the dimensions of charge, Equation [ VIII-27a ] is 



Since 



3 

usually written 

pJ'({fi*''})-pf+L4(^)<' + -- [VIII-27a'] 

where e*p(5)isthe dynamic effective charge tensor associated with the s 

nucleus. The nuclear displacements are a linear combination of the normal modes of 
the system 

Up'=E rp(5;v)^v) [Vni-28] 

v=l 

so that we may recast Equations [ VIII- 27a ] and [ VIII- 27b ] in the form 

pj' ({fi*1)-/^r + i:4(v) ^v) + - [ VIII-29a] 



v=l 

N 3N 

P?=L4;yEp ui^'=i;ai5'(v)Ep ^v) [VIII-29b] 

i=l v=l 



N N 

where el{v)=Y,el^{s) rp(^; v) and a^^^ (v)= J^a^^'^ r^(^;v) 

1=1 s=l 
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Quantum Theory of Raman Scattering 

The Raman nonlinearity provides a coupling between the photon field described by the 
Hamiltonian 



^ rod =11 ^(O-k 



Ha 



a-, a: +- 

k,a k, o 



[ VIII-30 ] 



and the phonon field described by the Hamiltonian 



[ vni-31 ] 



where b ^ and b ^ are, respectively, the phonon creation and destruction operators. In 
the electric dipole approximation, the interaction Hamiltonian can be written 



, 3^ 

^int=-2WW EaEp ^(V) 



v=l 



[ VIII-32 ] 



In earlier discussions we have shown that 



n 



[ VIII-33a ] 



and 



1/2 



V 



Mk,a)(a;^^-atJ [ VIII- 33b ] 



so that the complete interaction Hamiltonian may expressed as 
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271 n 



si; 2 



v=i V 



[ VIII-34 ] 



k', a' 



To concentrate on the key issue, we write scattering Hamiltonian which is 
proportional to o.^^^ a^,^, and which describes the scattering of a photon from the 

{k',a'} (initial) state to the {k,a } (final) state 



where a^^\l ,F;v)= al^{v) e^{l) eJP) 



2% n 

V 



n 



.1/2 



2M(0 



((0, CO^)'^' a^''\l,F;v) a, >^ [ VIII-35 ] 



V J 



We may now use the Fermi golden rule'*'' to calculate the transition rate per 

"molecule" for inelastic Stokes scattering events ~ viz. for 
{nj,np,n^} ^{n^-l,n^ + l,n^ + l} 



2 



n 



2% n 

V 



2M(0^. 



((0, la^^'Xl ,F;v)\ 



[ VIII-36a ] 



1 Stt' 



n 



2M(0 



V J 



(cD, C0^)| ,F; v) I' (n^ + 1) (n^ + 1) 5(a), -(0 ^ -(0 J [ VIII-36b ] 



See the appendix to this section for a recapitulation of the theory of transition rates. 
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and the transition rate per "molecule" for inelastic anti-Stokes scattering events 
viz. for [n,, np,i\]^{nj -I, Hp + l,n^ -1} 



1 27C 



87C' 



27t n 



2M (0 



V J 



(CO, co^)|fl<^'(/,F;v)|' 



x|(n,-l,n^ + l,n^-l| a'l a , B ^\n, ,np,n^f b{(£)j +(£)^-(£) p) [ VIII-36c ] 



2M(0 



V J 



{(£), CO^)|a*^'(/,F;v)|^ n, (n^ + 1) 5((0, +(0, -(0^) 



For a total number of "molecules" 1^p^„i we may write the total time rate of change of 
photons of frequency (0^ as 



d n. 



d t 



d Up 



87t P„,.,CD,a), 
V 



y2M a),y 



a'"'(/,F;v) 



Xti; {rip + l) (n^ + l)5(cO/ -COf-O) J 



d t 



AS 



V 



y2M CO,y 



XHj {np + l) n^5((0,+C0^-C0^) 



[ VIII-37a ] 



[ VIII-37b ] 



Stimulated Raman Effect: Simple model 

Suppose that a high power laser at (0, — the pump -- propagates parallel to the z-axis 
and that the Raman medium occupies the half-space z > . The spatial rate of change 
of the Stokes photons at co^ — the signal -- may be written 



50 



^ = G,n,{z)[np{z)+l] 

dz 



[ VIII-38 ] 



Assuming for the nonce that ~ . 
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where 



cV y. 



2M CO,, 



[ VIII-39 ] 



V y^'" J 



0.20 



40 W/cm^ 
60 W/cm^£ 
100 Wjfcm^S 
a)0 W.-tiTfS 




-10 -5 5 10 15 20 £5 30 

Probe Detuning (MHz) 

Raman Gain in rubidium vapor for different coupling laser intensities 



The rate constant is included this expression for gain to account for the broadening 
of a given vibrational state. Each Stokes event, of course, also depletes the incident 
beam so that 



— 7^ =-Gr nj{z) [np{z)+l] 
a z 



[ VIII-40 ] 
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which means that nj,[z) + nj{z) = n^(0) . [ Vni-41 ] 

Substituting this conservation condition into Equation [ VIII-SS ] we see that 

= [n,(0) - n,{z)] [n,{z) + l] [ VIII-42a ] 

a Z 



or that 



Therefore 



d~ nJz) dnJz) r / ^^ 

r t Ut - r M law = ^ + ^ ^ ' ^ "^^^'"^^^ ^ 



exp{G^[l + n,(0)]z}-l 

nJz) = nAO) !^ j r i ^ [ Vni-43 ] 

^ ^,(0) + exp{G,[l + n,(0)]z} 

This simple model predicts that Stokes photons initially increase Unearly with distance 
and that there complete conversion at distances large compared to | [l+ n^(0)] | 



I.e. 



lim np{z) =G^[l + nj{0)] z [ Vni-44a ] 

lim np{z) = n,{0) [ VIII-44b ] 
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n|(0) = 10 




0.4 O.S 1.2 1.6 



Coherent anti-Stokes Raman Spectroscopy (CARS): 

There is an important four wave mixing process involving the Raman nonlinearity. We 
may see how this comes about by realizing that Equation [ VIII-29 ] implies a potential 
energy contribution 

3N 

:;L<'(v)E„Ep^(v) [Vlll-45 

^ v=l 

so that inclusion of the Raman nonlinearity modifies the equation of motion for the 
vibrational normal modes - viz. 

q{t,v)+C0l q{t,v)= ^ a5'(v) E„ Ep [ Vni-46 ] 

When the normal mode is driven by a two (visible) wave input (see discussion of four- 
wave processes above), its dominant response will be at the difference frequency 
COJ-CO2 — viz. 
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^(^V) +(0l q{t,v) 

I [ vni-47 ] 

= — «ap (v) E„(co,) Ep(c02)exp[-/ {(0-(0,) t\ + c.c. 



so that 



. X 1 «1r'(v)£„(co,)£p*(co2) exp[-/(c0i-C02) 
q{t,v) = ' V ; — 2 — [ VIII-48 ] 

From Equation [ VIII-29b ] we see that excitation of a vibrational mode, through the 
Raman nonlinearity, generates a nonlinear polarization 



1 V ''V(^) «V(v)Ep((0,) E^((0,) E,{(0,) exp[-i [M-CD,) t] 
ra= — L —r ^2 [ VIII-49 ] 



3N 



at the four wave mixing frequency 20^-0) 2- Thus, we see that the medium has an 
effective third order susceptibility 

which exhibits a resonance when CO; -CO 2 equals (0^. 
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VIIIA. Appendix: Radiative Transition Rates Revisited 

This appendix builds on the formulation presented in Review of Basic Quantum 

Mechanics: Dynamic Behavior of Quantum Systems, Section II of the lecture notes 

entitled The Interaction of Radiation and Matter: Semiclassical Theory (hereafter 

referred to as IRM:ST) to obtain expUcit and reasonably general expressions for 

radiative transition rates. Let us suppose that a coupled system of radiation and matter 
is described in the Schrodinger picture by a wave function |^^(?)). According to 

Equation [ 11-23 ] of IRM:ST, we may write 




where ^hC^^ is the Hamiltonian of the uncoupled radiation and matter systems. Further, 
Equations [ I-36a ] and [ I-37a ] of IRM:ST inform us that 




[ VIIIA-2 ] 



where 




[ VIIIA-3 ] 
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M%{t)f = \M'ro{t)\%{t))f 

= l^Vf\'ro{t)'T,{tj,)<T,-\t,)\^,{t,)J 



[ VIIIA-4 ] 



To be more specific, if the system happened to be in state |\|/,) at a time t^, the 
probability that the coupled system is in a state |\|/^) at a time t is given by 



[ VIIIA-5 ] 



In the most experimental circumstances, we are interested in the rate at which 
transitions take place from some particular initial state \\\f.) to a set of final states | Vy) - 



I.e. 



I^tX (^'^o) = tX I\V/Ko(0 "^iM To-^o) |¥,)f [ VIIIA-6 ] 



T dt 



dt 



The first order approximation ( for 1/t ): 

Using the second term in the iteration set forth in Equation [ TL-llc ] of IRM:ST, we can write 



f^M^'O = !i2}UKo(0 |[-0 df<T-\t) exp(8 t) ^o(0|^oiO) 
■ Jexp(-i (0^ t) (v|/^ I v|/,.) j dt' exp{[t ((o^ - w ,) +8] r j 



= Km 



= Um 

£^0 



[ VIIIA-7 ] 



^ mt exp(e r-to), t) 

(o), -COy + i£) 



The Interaction of Radiation and Matter: Quantum Theory 



where the factor exp(8 t) was introduced to avoid the transient effects which might 

otherwise result from an apparent sudden application of the interaction between 
systems. Therefore, the first order approximation for l/x is given by 



ex 



p(2er) 



+e 



2k 



[ VIIIA-8 ] 



= T2 L|(¥/Kintl¥,)| S(o>,-o)^) 



which is, of course, the generalized form of the Fermi golden rule. 
The second order approximation ( for l/x ): 

Using the third term in the iteration set forth in Equation [ I-27c ] of IRM:ST, we can 
write 



£->0 



exp(E^^£M 



+ (¥/|'ro(0|[-^) Iclt'cxi^et) [VIIIA-9a] 
x}jrexp(8 0^H„t(0 ko"'(0) |x|/,> 



or in the Schrodinger picture 
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E->0 



+ (¥/|'ro(r)|[-^] l^^'MOexp(e0^int^o(0 [VmA-9b] 
X } rff T„-'(0 exp(E ^ "TolO ko"(0) I ¥.) 



Using the closure theorem we see that 



+ 



^ X exp(-t 0), (v, I V|/ J <X|/„ I |vi/,.) [ VIIIA-10 ] 



X I i/r' j dt" exp[i (0^ r'+E f- i (£>„{t'- t") +e t"- i CO,. 



Integrating and writing the matrix elements in more concise notation, we see that 



exp(E.-fco,..) [fn 



a),-a)„ + t£/2 



[VIIIA-11] 



Therefore, following the arguments presented above, we see that the second order 
approximation for l/x is given by 
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_2Tly 



to -to 



[ VIIIA-12 ] 
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The nth order approximation ( for l/x ): 

We then may make the obvious extrapolation of these results and write the nth order 
approximation for 1/ x as 



1 _ 271 y 



,f\nr |A,lv (/Kin.h)(^KintlO , 



/ 



+ 



ff^ m2 m^_j 



((0,-(0„ )((0,-(0„ )...((0,-(0„ ) 



[ VIIIA-13 ] 



5 (to, -to,) 



where \in^/, \in^/, etc. could be real imtermediate states or so called virtual 
intermediate states for the transitions to the state |/). The conservation of 
energy holds only between |/) and |/). 



